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SUMMARY 

l 

The  transonic  equivalence  rule  with  lift  has  been  studied  experimentally. 
Under  the  conditions  defined  by  the  rule,  wing-body  models  with  the  same 
cross-sectional  area  distribution  but  with  wings  of  different  aspect  ratios  and 
thicknesses  were  tested  at  transonic  speeds.  The  correlations  of  the  lift 
depending  wave  drag  and  the  lift  parameter  for  the  models  are  found  to  be 
similar.  The  similitude  of  the  outer  flow  field  under  the  same  lifting  condi¬ 
tion  is  also  established.  The  drag-rise  of  the  equivalent  body  of  revolution 
with  the  cross-sectional  area  distribution  including  that  due  to  lift  correlates 
well  with  those  of  the  wing-body  models  at  the  design  lift  condition.  The 
experiment  therefore  verifies  the  transonic  equivalence  rule  at  lifting  condi¬ 
tion. 


RESUME 


On  a  etudie  experimentalement  l’application  a  la  portance  de  la  regie 
d’equivalence  transsonique.  Dans  les  conditions  definies  par  la  regie,  des 
profils  aerodynamiques  ayant  la  meme  distribution  de  section  transversale 
mais  des  allongements  et  des  epaisseurs  differents  ont  ete  soumis  a  des 
essais  a  des  vitesses  transsoniques.  On  a  trouve  pour  le  parametre  de  portance 
et  la  trainee  d’onde,  liee  a  la  portance,  des  correlations  similaires  entre  les 
modeles  consideres.  On  a  aussi  etabli  la  similitude  des  champs  d’ecoulement 
exterieurs  dans  les  memes  conditions  de  portance.  II  y  a  une  bonne  corre¬ 
lation  entre  l’augmentation  de  trainee  du  corps  de  revolution  equivalent,  en 
ce  qui  a  trait  a  la  distribution  de  section  transversale  incluant  celle  due  a  la 
portance,  et  celle  des  profils,  aux  conditions  de  portance  theorique.  L’expe- 
rience  permet  done  de  verifier  la  validite  de  la  regie  d’equivalence  transso¬ 
nique  aux  conditions  de  portance. 
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AN  EXPERIMENTAL  STUDY  OF 
THE  TRANSONIC  EQUIVALENCE  RULE  WITH  LIFT 


1.0  INTRODUCTION 

In  the  last  two  decades,  the  transonic  area  rule  has  played  an  important  role  in  the  design  of 
high  speed  aircraft.  The  rule  indicates  the  relation  between  the  transonic  drag  rise  and  the  cross- 
sectional  area  distribution  of  the  aircraft,  allowing  the  configuration  to  be  designed  with  niminal  drag 
rise.  The  principle  of  the  rule  is  based  on  the  properties  of  the  air  moving  at  a  speed  near  the  speed  of 
sound.  At  Mach  number  close  to  unity,  the  air  flow  will  resist  any  streamtube  area  change  with  a  large 
force.  To  avoid  large  perturbation  pressures  and  hence  high  drag,  the  changes  of  the  cross-sectional 
area  should  therefore  be  gradual  and  smooth.  The  constraint  of  the  stream-tube  effectively  causes  the 
cross-flow  in  planes  normal  to  the  free-stream  direction  to  behave  like  an  imcompressible  flow.  This 
allows  the  nonlinear  transonic  flow  field  to  be  analyzed  asymptotically  leading  to  the  formulation  of 
the  area  rule. 

The  analysis  shows  that  the  far  field  is  governed  by  the  transonic  small  disturbance  theory. 
The  cross-flow  near  the  body  tends  to  be  incompressible  as  in  the  case  of  the  flow  near  a  slender  body 
and  hence  can  be  treated  in  a  similar  manner.  To  the  far  field,  the  displacement  effect  of  the  body 
appears  as  that  generated  by  a  line  source  distributed  along  the  axis  of  the  body  and  is  thus  identical 
to  that  of  an  axisymmetrical  body  having  the  same  cross-sectional  area  distribution.  This  similitude  of 
the  displacement  effect  is  basically  the  equivalence  rule  or  the  area  rule*1,2) .  The  concept  was  verified 
experimentally  by  Whitcomb13  *  and  has  since  been  applied  to  aircraft  design  for  drag  reduction  at 
transonic  speeds. 

When  the  aircraft  is  in  lifting  condition  the  transonic  flow  past  the  wing  causes  further 
expansion  of  the  streamtube.  To  the  far  field  the  lift  distribution  along  the  body  axis  is  represented 
by  a  line  doublet  and  the  expansion  of  the  streamtube  appears  as  additional  sources  which  relate 
directly  to  the  lift  distribution.  To  form  the  equivalence  body,  the  effective  area  representing  these 
additional  sources  must  be  added  to  the  geometric  cross-sectional  area.  The  basic  idea  was  incorpo¬ 
rated  in  the  design  of  an  early  experiment  of  transonic  transport  aircraft,  but  the  amount  of  lift  com¬ 
pensation  area  was  then  derived  entirely  empirically.  Theoretical  analysis  of  this  problem  was  later 
carried  out  by  Cheng  and  Hafez*5*,  and  Barnwell16,7*.  In  their  analysis  the  detailed  structure  of  the 
flow  is  delineated  and  the  conditions  of  the  similitude  of  the  far  fields  are  established.  The  explicit 
expressions  of  the  doublet  and  source  strengths  induced  by  lift  given  in  the  analysis  are  particularly 
important  for  engineering  applications.  This  has  been  demonstrated  in  a  design  study  of  transonic 
transport  aircraft*8*,  and  the  estimations  of  drag-rise  for  configurations  of  fighter-type  aircraft*9*. 

The  extended  equivalence  rule  with  lift  can  be  expressed  in  a  form  of  flow  similitude.  It 
states  that  flows  having  the  same  distributions  of  sources  S',(x)  and  doublets  D0(x),  with  the  same 
transonic  similarity  parameter  K  are  equivalent.  The  nonlinear  structure  of  the  flow  field  including 
the  formation  of  shock-waves  and  the  associated  drag-rise  will  be  identical.  The  parameters  controlling 
the  lift  effect  and  the  geometry  of  the  body  do  not  appear  explicitly  and  it  is  therefore  more  general 
than  the  classical  transonic  similarity  law*10*. 

The  present  study  attempts  to  establish  the  similitude  of  the  flow  fields  according  to  the 
extended  equivalence  rule.  Two  wing-body  models  are  designed  with  the  wings  having  different 
thicknesses  and  aspect  ratios  but  the  same  cross-sectional  area  distributions.  The  models  are  tested 
at  incidences  through  the  transonic  Mach  number  range.  The  lift  and  the  drag  forces  acting  on  the 
models  are  measured  and  the  static  pressure  variation  at  the  far  field  is  also  monitored  in  the  experi¬ 
ment.  The  data  are  then  examined  in  the  parametric  forms  defined  by  the  equivalence  rule.  The  main 
parameter  chosen  for  detailed  study  is  the  lift-depending  drag  associated  with  the  shock  wave  forma¬ 
tion.  The  result  shows  that  for  the  same  transonic  condition,  the  correlations  of  the  wave  drag  and  the 
lift  parameter  for  the  two  models  are  similar.  Comparisons  of  far  field  pressure  distributions  and  shock 
patterns  for  the  two  models  also  indicate  similarity  of  the  flow  fields.  Thus  the  flow  similitude  defined 
by  the  equivalence  rule  is  obtained. 


Two  equivalent  bodies  of  revolution  corresponding  to  the  conditions  of  zero  lift  and  of 
finite  lift  respectively  are  also  part  of  the  experiment.  In  both  cases  the  drag-rise  of  the  equivalent 
bodies  correlates  very  well  with  those  of  the  wing-bodies  at  the  design  conditions.  The  latter  result 
verifies  the  effective  displacement  of  the  flow  due  to  lift  as  defined  by  the  extended  equivalence  rule. 

In  the  report,  the  theoretical  results  of  the  extended  equivalence  rule  is  first  outlined.  The 
design  of  the  models  for  the  experiment  is  then  described.  The  experimental  procedure  and  the  data 
are  presented  in  some  details.  Finally,  the  data  are  analyzed  in  the  parametric  forms  defined  by  the 
equivalence  rule  and  the  conditions  of  flow  similitude  are  then  established. 


2.0  TRANSONIC  EQUIVALENCE  RULE  WITH  LIFT 

Transonic  flows  past  a  three-dimensional  configuration  have  been  analyzed  by  Cheng  and 
Hafez,  leading  to  the  establishment  of  the  equivalence  rule  with  lift151.  The  analysis  shows  that  at  a 
specified  transonic  similarity  parameter  K,  the  outer  flow  is  uniquely  determined  by  a  line  doublet 
distribution  D0(x)  due  to  lift  and  a  line  source  distribution  S^,(x)  related  to  both  lift  and  displace¬ 
ment  of  the  body.  Flows  having  the  same  distributions  of  D0(x)  and  Se(x),  with  the  same  value  of  K 
are  therefore  equivalent.  These  parameters  are  given  as  follows: 


Mi  -  1 
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The  basic  parameters  are  defined  as 
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and  the  outer  region  of  the  flow  is  scaled  as 
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where  the  superscript  *  denotes  quantities  with  physical  dimensions,  V.  is  the  body  length,  b  the 
semispan  of  the  wing  and  t  the  ratio  of  the  transverse  length  scales  for  the  inner  and  the  outer  regions 
of  the  flow  field.  S*m  and  F  *m  are  the  maximum  cross-sectional  area  and  the  maximum  of  the 

local  lift  force  F  (x)  respectively.  Thus  r  is  the  parameter  of  thickness  of  the  configuration,  a,  the 
lift  and  X  the  wing  sweep  parameters.  The  parameters  controlling  the  lift  effects  on  the  outer  flow  are 
v.  and  l\  .  D()(x)  is  the  line  doublet  due  to  lift  and  Sc(x)  is  the  line  source  in  the  form  of  the  first 
derivative  of  the  effective  cross-sectional  area  consisting  of  the  geometric  cross-sectional  area  S  is) 
and  the  contributing  due  to  lift.  The  first  term  of  the  lift  effect  in  the  square  bracket  of  Equation  (3) 
is  from  the  axial  distribution  of  the  lift  force  F(x).  The  last  term  E(x)  signifies  the  cross-flow  kinetic 
energy,  comparable  to  induced  drag  and  is  finite  in  the  wake  region. 

When  the  flow  fields  are  equivalent,  the  local  pressure  coefficient  and  the  Mach  number 
are  functions  of  the  transonic  similarity  parameter  K  only 
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rX 


f  (x,  rj ,  co;  K) 


Hence  the  drag-rise  associated  with  the  shock  waves  is  also  a  function  of  K 


(10) 


Dtt/(P„Uib-  r2  M£)  =  f  (K)  (11) 

The  conditions  of  the  equivalence  rule  as  given  in  Equations  (1),  (2)  and  (3)  are  much  less 
restricted  than  the  classical  transonic  similitude.  The  parameters  controlling  the  lift  o,  and  the  wing 
geometries  for  the  equivalent  flows  do  not  require  to  be  the  same. 


We  have  summarized  the  results  of  the  extended  equivalence  rule  with  lift  as  given  in 
Reference  5.  The  parameters  shown  in  Equations  (1)  to  (4)  are  complicated  in  form,  but  it  is  not 
difficult  to  trace  their  relations  to  the  respective  physical  quantities.  The  basic  form  of  these  param¬ 
eters  can  be  derived  from  the  similarity  analysis  of  the  transonic  small  disturbance  equation.  This 
exercise  does  help  us  to  understand  the  formation  of  these  parameters  and  their  physical  implications 
and  is  therefore  included  in  the  report  as  Appendix  A. 


3.0  DESIGN  OF  THE  EXPERIMENT 

The  purpose  of  the  experiment  is  to  establish  the  flow  field  similarity  according  to  the 
equivalence  rule  at  lifting  condition.  This  is  done  by  correlating  the  flow  parameters  in  the  form 
defined  by  the  equivalence  rule  and  if  the  correlations  for  different  models  are  similar,  then  the  flow 
fields  are  equivalent.  With  the  understanding  that  the  ultimate  application  of  the  equivalence  rule  is 
to  reduce  the  drag-rise  caused  by  the  formation  of  shock-waves,  the  lift  dependent  wave  drag  is 
chosen  as  the  main  parameter  for  the  present  study.  Its  direct  association  with  the  shock  formation 
in  the  outer  field  provides  a  measure  of  the  nonlinear  features  in  the  region.  In  addition  the  static 
pressure  variation  in  the  far  field  is  also  monitored  to  add  further  information  for  the  establishment 
of  flow  similitude. 

Two  wing-body  models  are  designed  for  the  experiment.  The  wings  of  the  models  have 
different  thickness  and  aspect  ratio  and  are  mounted  at  the  central  plane  of  a  cylindrical  body  housing 
an  internal  balance  for  force  measurements.  In  the  design  of  the  models  the  functual  forms  relating 
to  the  parameters  given  in  Equations  (2)  and  (3)  must  be  considered.  These  relations  are  now  written 
in  terms  of  the  similarity  parameters  discussed  in  the  Appendix  A  for  a  better  understanding  of  their 
physical  implications. 

«s 

Do(x)  =  e  —  F(x) 


(12) 
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Sc(x)  =  St(x) 
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The  sectional  profile  of  the  wings  is  chosen  symmetrical  so  that  the  lift  effect  can  be  simply 
related  to  the  angle  of  attack.  The  wing  planforms  are  selected  to  be  geometrically  similar  so  that  the 
functions  F(x),  T(x)  and  E(x)  are  the  same  for  both  models.  If  the  planforms  are  not  geometrically 
similar,  it  is  necessary  to  solve  inversely  from  Equations  (5)  to  (7)  for  the  planform  shape  of  the 
second  model  if  the  first  one  is  initially  specified.  For  geometrically  similar  planforms  the  condition 
D0(x)  being  the  same  is  satisfied,  if  the  parameter  e  ajT  is  identical  for  both  models  (see  Eq.  (12)). 
The  condition  Se(x)  being  identical  requires  that  the  parameters  e  and  X  are  the  same  for  the  two 
models  respectively.  This  leads  to  the  trivial  case  that  the  two  planforms  are  exactly  identical.  How¬ 
ever,  if  the  function  Fx  is  proportional  to  the  function  E(x),  then  a  functional  relationship  can  be 
established  for  e  and  X  satisfying  the  desired  similarity  condition  of  Se(x).  The  effective  areas  can  thus 
be  kept  identical  for  the  two  models  but  with  different  sweep  back  parameters.  From  the  slender  wing 
theory  it  can  be  shown  that  only  a  delta  planform  satisfied  this  condition.  Thus  the  delta  planform  is 
adopted  for  the  model  design.  The  evaluations  of  the  functions  F(x),  T(x)  and  E(x)  for  a  planar  wing 
at  angles  of  attack  are  given  in  Appendix  B  and  the  results  are  summarized  as  follows: 


F(x) 

Fx2(x) 

T(x) 

E(x) 


4  a"  aj 


(14) 


where  a(x)  is  the  leading  edge  contour.  The  expressions  have  been  normalized  by  the  value  of  F(x)  at 
the  trailing  edge.  For  a  delta  wing  ax  is  unity  and  Fjj  is  thus  proportional  to  E(x). 

3.1  Planform  Design 

The  respective  relations  of  X,  r,  a  and  for  two  wings  are  now  determined  in  the  following 
procedure.  The  geometric  cross-section  Sc(x)  is  made  identical  for  both  models,  since  the  lengths  of 
the  models  are  the  same.  Thus  from  the  definitions  in  Equation  (8),  we  have 

T,x,  =  r2x2  (15) 

which  determines  the  relative  thickness  and  aspect  ratio  of  the  wings.  The  condition  that  the  transonic 
similarity  parameter  K  being  identical  leads  to  the  Mach  number  relation 

M  = 

~t 

From  the  definition  of  e,  we  thus  have 

CJ_  _ 

f 2  ^2 


M 
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The  requirement  that  the  doubled  strength  distribution  Do(x)  is  the  same  implies  from  Equation  (12) 

that  e  —  assumes  the  same  value  for  both  wings, 
r 


(18) 
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With  this  condition  and  the  properties  of  delta  wing  given  in  Equation  (14)  the  requirement  that  the 
effective  area  induced  by  lift  is  the  same  leads,  from  Equation  (13),  to  the  relation 


With  Equation  (17)  a  relation  between  X,  and  X2  is  obtained 


(20) 


In  the  design  Xi  is  specified  and  X2  is  determined  from  Equation  (20).  The  relation  between  r,  and  r2 
follows  from  Equation  (15)  and  between  a6l  and  ^  from  Equation  (18). 


3.2  Wing  Section 

A  symmetric  airfoil  section  based  on  Boerstoel’s  designs*1 3'14)  is  used  for  the  wing  section. 
The  nose  shape  of  the  original  section  has  been  modified  to  eliminate  the  peaky  pressure  rise*15*.  The 
sectional  shape  is  scaled  to  the  required  thickness  of  the  model.  The  coordinates  of  the  airfoil  is  given 
in  Table  1  for  thickness-chord  ratio  of  0.1.  Typical  calculated  pressure  distribution  on  the  upper  and 
the  lower  surface  of  the  airfoil  at  angle  of  attack  is  shown  in  Figure  1.  The  pressure  near  the  nose 
region  decreases  monotonically  and  there  is  no  pronounced  suction  peak.  This  condition  is  required 
for  the  existence  of  the  solution,  which  leads  to  the  equivalence  rule*5'8). 


3.3  Wing  Tip  Round  Off 

The  equivalent  body  shape  as  given  in  Equation  (13)  depends  on  the  functions  Fx ,  E(x)  and 
T(x)  when  the  model  is  in  lifting  condition.  For  a  delta  wing,  all  these  functions  will  be  discontinuous 
at  the  trailing  edge  of  the  wing  leading  to  a  discontinuity  of  the  equivalent  body  contour.  In  order  to 
smooth  the  equivalent  body  shape  (as  the  trailing  edge  position  is  approached)  to  prevent  adverse 
effects  on  the  boundary  layer  flow,  the  wing  tip  is  shaped  as  a  sine  function  to  allow  the  slope  of  the 
leading  edge  contour  to  decrease  gradually  to  zero .  The  contour  is  given  as 


0.9  <-  <  1, 
c 

0.9  <  y/b 


(21) 


The  contouring  of  the  wing  tip  is  in  contradiction  to  the  similarity  requirement  as  discussed  at  the 
beginning  of  Section  3.  However,  since  only  a  small  portion  of  the  wing  planform  has  been  altered, 
the  effect  to  the  outer  flow  field  is  small  and  the  similarity  of  the  flow  can  still  be  established. 


3.4  Wing-body  Combination 

For  a  wing-body  combination  the  parameter  scaled  to  the  wing-body  can  be  related  to  those 
scaled  to  the  wing  alone*8*.  If  the  length  scale  2  denotes  the  body  length  and  rQ  is  the  root  chord  of 
the  wing,  thus  the  following  relations  can  be  established 
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and  with 


x  -  x 
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/  £  \2  _ 
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Equations  (12)  and  (13)  become 
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(24) 
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where  x0  is  the  location  of  the  apex  of  the  wing  planform  along  the  body  axis.  The  expression  inside 
the  square  bracket  of  Equation  (25)  is  the  same  as  the  wing  alone  case  Equation  (13),  except  that  the 
scaling  parameter  is  now  el  The  conditions  discusses  in  Section  3.1  can  be  applied  to  the  wing-body 
case  as 


X,  (26) 


and  the  X  relation  in  Equation  (20)  still  applies.  The  wings  can  therefore  be  adopted  to  the  wing-body 
model  without  alteration. 


The  nose  contour  of  the  central  body  is  a  von  Karman  ogive  which  gives  minimum  drag  for 
a  given  base  area  at  supersonic  speed.  The  nose  shape  is  not  optimized  at  transonic  speeds  and  a  shock 
will  be  formed  at  the  end  of  the  nose  portion. 


The  dimensions  and  the  parameters  of  two  models  are  given  in  Table  2.  The  general  layout 
of  the  models  is  shown  in  Figure  2.  Since  model  WB1  has  a  thicker  wing  and  a  larger  portion  of  the 
wing  relative  to  the  span  is  covered  by  the  central  body  the  cross-section  area  of  the  wing  (excluding 
the  part  inside  the  body)  is  smaller  than  that  of  model  WB2  at  the  corresponding  section.  The  central 
body  of  the  model  WB1  is  thus  enlarged  slightly  to  compensate  for  the  area  defect  of  the  wing  so  that 
the  total  area  of  the  wing-body  is  the  same  for  both  models.  The  cross-sectional  area  distribution  is 
shown  in  Figure  3. 


For  a  configuration  with  thin  wing  and  slender  body  the  classical  area  rule  is  applicable  for 

/  tV/3 

the  parameter  2X1 —J  less  than  unity'  '*  n .  The  present  models  have  values  of  0.2  and  0.3  respec¬ 
tively  and  are  well  within  the  range  of  application. 


3.5  Equivalent  Bodies 


At  zero  lift,  an  equivalent  body  can  be  simply  designed  to  have  the  same  cross-sectional 
area  distribution  as  the  wing-body  models.  An  axisymmetric  model,  Bl,  has  been  designed  for  this 
condition  to  check  whether  the  classical  area  rule  applies  for  the  present  models.  When  the  wing-body 
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is  at  lifting  condition,  the  effective  area  due  to  lift  as  given  by  the  right-hand  side  of  Equation  (25) 
should  be  added  to  the  geometric  area.  The  parameter  controlling  the  lift  is  <f3j,/?'as  indicated  in 
Equation  (24).  An  equivalent  body  can  therefore  be  designed  for  a  lifting  condition  with  specific 
«Tt\  It' and  el  To  study  the  drag-rise  at  a  specified  lift  condition,  the  lift  parameterTcq  /? must  be  kept 
constant  as  Mach  number  increases.  However,  'e  is  a  function  of  Mach  number  and  the  effective  area 
defined  in  Equation  (25)  depends  one  Thus  a  series  of  axisymmetric  equivalent  bodies  would  be 
needed  to  cover  a  Mach  number  range  for  the  lifting  condition.  We  also  note  that  c~ appears  in  the  form 
|  Kn  e  |  in  Equation  (25)  and  within  the  range  of  Mach  number  of  0.8  to  1 .1,  i  fcn  e*  |  varies  less  than 
10  percent.  The  corresponding  effect  on  the  maximum  cross-sectional  area  is  only  2.5  percent.  Thus 
for  this  experiment  only  one  axisymmetrical  body  is  designed  at  a  specified  lift  condition  and  the 
same  model  is  used  to  study  the  drag-rise  through  the  Mach  number  range.  The  conditions  chosen  for 
the  design  are  as  follows. 


eat 


0.2024, 


e,  =  0.0381 


which  corresponds  to  the  wing-body  model  WB1  at  the  conditions 

M  =  0.975,  C.  =  0.37 

oo  L 

These  conditions  are  chosen  to  simulate  a  lifting  condition  at  cruise  of  a  typical  transonic  aircraft 
(see  Ref.  8).  The  angle  of  attack  of  the  aircraft  would  be  5  degrees  based  on  slender  wing  theory.  The 
effective  area  due  to  lift  is  also  shown  in  Figure  2.  The  added  area  on  the  aft-body  is  caused  by  the 
cross-flow  kinetic  energy  in  the  wake,  a  finite  contribution  from  the  function  E(x). 

The  models  are  designed  to  have  a  central  basic  body  for  housing  the  internal  balance.  An 
outer  body  consisting  of  the  wing-body  or  the  equivalent  body  is  slid  over  the  central  body  to  provide 
the  proper  contour.  For  the  lifting  equivalent  body  a  sleeve  is  slid  over  the  aft-body  to  increase  the 
cross-sectional  area  as  required.  A  nose  portion  common  to  all  models  is  then  installed  to  complete  the 
configuration.  The  general  layout  of  the  model  parts  before  assembly  is  shown  in  Figure  4. 


It  has  been  found  in  Reference  18  that  there  is  no  essential  difference  on  transonic  drag-rise 
of  the  wing  with  and  without  boundary  layer  tripping.  Thus  in  this  experiment,  performed  at  reason¬ 
ably  high  Reynolds  number,  no  boundary  layer  tripping  is  applied  on  the  wing  and  the  nose  of  the 
body  of  the  models. 

3.6  Flow  Similarity  Within  a  Wind  Tunnel 

Since  the  experiment  will  be  performed  in  a  wind  tunnel,  the  effect  of  the  tunnel  walls  to 
the  flow  past  the  model  must  be  considered.  The  severity  of  the  constraint  of  the  tunnel  wall  can  be 
estimated  by  comparing  the  physical  dimensions  of  the  working  section  and  the  length  scales  of  the 
flow  field.  The  transverse  length  scale  of  the  outer  region  of  the  flow  is  b/e" and  has  a  value  of 
106  inches  (2.69  m)  for  the  present  wing-body  model  and  the  design  conditions  discussed  in  the  last 
section.  The  wind  tunnel  has  a  60-inch  (1.5  m)  square  working  section.  Thus  the  tunnel  wall  inter¬ 
ference  cannot  be  treated  as  a  weak  correction  and  the  walls  should  be  considered  as  an  integral  part 
of  the  flow  boundary. 

For  the  flows  within  the  tunnel  to  be  similar  for  both  models  the  length  scales  in  both  the 
longitudinal  and  the  transverse  directions  must  be  the  same  respectively.  The  first  condition  is  satis¬ 
fied  since  the  models  have  the  same  length.  The  second  condition  follows  from  the  second  design 
condition  specified  in  Equation  (26),  giving  identical  transversal  length  scales.  Since  the  walls  are  now 
an  integral  boundary  of  the  flow,  the  boundary  condition  at  the  wall  must  satisfy  the  similarity 
requirement  as  well.  For  a  tunnel  with  perforated  walls  the  similarity  scaling  shows  that  the  porosity 
of  the  wall  should  be  the  same  for  the  tests  of  two  models.  This  condition  is  automatically  satisfied 
for  the  present  tunnel  with  fixed  wall  porosity.  A  similarity  scaling  of  the  boundary  condition  at  the 
wall  is  given  in  Appendix  C. 
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4.0  EXPERIMENTAL  PROGRAMS 

4.1  Experimental  Methods 

The  experiments  were  performed  in  the  NAE  5  ft  (1.5  m)  X  5  ft  (1.5  m)  blowndown  wind 
tunnel.  The  tunnel  system  and  performance  are  described  in  detail  in  Reference  20.  For  the  present 
experiment,  the  5-ft  (1.5  m)  square  transonic  section  with  perforated  walls  was  used.  The  model  was 
mounted  on  the  supporting  sting  which  could  be  pitched  at  the  desired  angles  with  the  model  remain¬ 
ing  in  the  center  of  the  Schlieren  window  at  the  transonic  downstream  position.  The  high  speed 
Sehlieren  system  was  in  operation  during  the  tests.  A  1-inch  (25.4  mm)  Task  Corporation  MK  XIV 
6-component  balance  was  used  for  measurements  of  forces  and  moments.  The  balance  was  pre¬ 
calibrated  outside  the  tunnel  before  the  test.  The  base  pressure  of  the  model  was  measured  by  a  pres¬ 
sure  transducer  installed  in  the  cavity  of  the  sting.  Static  pressure  pipes  of  1  inch  (25.4  mm)  diameter 
were  mounted  along  the  center  line  of  the  top  and  the  bottom  walls  respectively  for  measurements  of 
static  pressure  distribution  near  the  wall. 

Four  models  were  tested  in  the  program:  two  wing-bodies,  WB1  with  small  wing  and  WB2 
with  large  wing;  and  two  equivalent  bodies  of  revolution,  B1  for  zero  lift  and  B2  with  lift.  The  test 
Mach  numbers  ranged  from  0.8  to  1.1.  The  Reynolds  number  in  the  experiment  was  fixed  at  4  X  106 
based  on  the  root  chord  of  6  inches  (152.4  mm).  High  Reynolds  number  tests  were  also  performed  at 
8  X  106  for  Mach  number  of  0.94  for  some  models  to  check  the  Reynolds  number  effect.  During 
each  run,  the  wing-body  model  WB1  was  pitched  continuously  from  -4°  to  10°  at  a  pitch  rate 
1.887  deg/sec  and  then  followed  by  three  stepping  incidences  at  0°,  5°  and  10°  nominally.  During 
the  constant  incidence  steps,  enough  time  was  allowed  for  the  static  pressure  at  the  upper  and  lower 
walls  to  be  scanned.  Schlieren  photographs  of  the  flow  field  were  taken  during  the  test  at  a  rate  of 
0.8  sec/frame.  For  the  WB2  model  with  a  larger  wing  span,  the  range  of  incidence  was  half  of  that  of 
WB1  and  the  pitch  rate  was  also  reduced  accordingly. 

Tunnel  flow  angularity  was  checked  at  M„  =  0.8  by  testing  the  model  B1  in  the  normal  and 
the  inverted  positions.  The  flow  angularity  and  the  geometric  asymmetry  of  the  model  were  deduced 
from  these  tests.  The  geometric  asymmetry  of  the  model  and  the  flow  angularity  obtained  were  read 
in  directly  to  the  data  reduction  program  for  the  correction  of  the  geometric  incidence  of  the  model. 
The  sting  used  in  the  experiment,  unfortunately,  could  not  be  rotated  automatically  by  the  roll 
mechanism  of  the  supporting  system  and  the  model  had  to  be  inverted  manually.  This  was  a  time- 
consuming  process  and  therefore  the  flow  condition  was  checked  for  one  Mach  number,  =  0.8 
only.  The  flow  angularity  was  found  to  be  0.53°  and  the  asymmetry  of  the  model  caused  a  further 
deviation  of  0.25°.  Since  none  of  the  models  had  been  checked  for  inaccuracy  of  machining  due  to 
lack  of  time,  the  degree  of  asymmetry  might  be  different  between  the  models.  Thus  only  the  flow 
angle  offset  was  incorporated  in  the  data  reduction  program  and  the  value  of  0.55°  was  used  for  all 
models  and  Mach  numbers  in  the  tests. 

Because  of  uncertainty  of  the  flow  angularity  at  this  range  of  Mach  numbers,  the  axi- 
symmetrical  bodies  were  also  tested  in  a  similar  “ramp  and  step”  program  with  the  pitch  angle  ranged 
from  -3°  to  3°  and  steps  at  0.1°,  0.2°  and  0.3°,  from  which  the  minimum  drag  was  obtained. 

Surface  flow  visualization  was  applied  to  a  few  tests  to  provide  information  on  boundary 
layer  transition,  separation  and  formation  of  leading  edge  vortices. 

4.2  Data  Reduction 

The  experimental  data  were  reduced  on-line  after  each  run.  Some  important  data  such  as 
CL  vs  ck  and  CD  vs  C,  were  plotted  after  the  data  had  been  reduced  to  coefficient  forms.  This  large 
amount  of  data  were  further  analyzed  using  the  computer  of  the  NRC  Computation  Center  and 
organized  in  parametric  forms  suitable  for  the  present  studies. 

The  preliminary  results  of  the  wing-body  models  indicate  that  the  CL  vs  a  curves  do  not 
pass  through  the  origin  of  the  co-ordinates.  The  axis  of  symmetry  of  the  CD  vs  CL  polars  also  offset 
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slightly  from  the  CL  axis.  This,  as  discussed  in  the  previous  section,  is  due  to  the  asymmetry  of  the 
model  and  that  the  tunnel  flow  angularity  at  other  Mach  numbers  may  be  slightly  different  for  that  of 
Mx  =  0.8  which  is  used  in  the  data  reduction  program  for  flow  angularity  correction.  These  effects 
are  corrected  for  the  experimental  data  and  both  CL  vs  a  and  Cu  vs  CL  curves  are  adjusted  to  appear 
as  from  a  symmetric  model  lining  up  properly  with  the  tunnel  flow.  The  estimation  of  the  offsetting 
of  the  data  is  given  in  Appendix  D.  The  data  are  then  presented  in  terms  of  CL ,  a  and  Cu  in  Fig¬ 
ures  5  and  6  for  the  models  WB1  and  WB2  respectively.  The  values  of  these  data  are  listed  in  Tables  1 
and  2  respectively. 

The  data  presented  have  not  been  corrected  for  the  wind  tunnel  wall  interference.  It  has 
been  discussed  in  Section  3.6  that  at  transonic  speeds,  the  tunnel  wall  interference  can  no  longer  be 
treated  as  a  weak  correction.  Thus  the  data  are  presented  as  obtained  with  the  tunnel  walls  as  an 
integral  part  of  the  flow  boundaries.  It  is  desirable,  nevertheless,  if  the  effect  of  the  wall  constraint 
can  be  estimated  even  in  the  low  Mach  number  range.  In  the  limit,  the  subsonic  theory  may  be  con¬ 
sidered  to  be  applicable  asymptotically.  The  subsonic  interference  theory  shows  that  the  angle  of 
attack  and  the  drag  corrections  are  proportional  to  the  area  ratio  of  the  model’s  maximum  cross- 
section  to  that  of  the  working  section  of  the  tunnel*2 **•  For  the  present  experiment,  the  ratio  has  a 
value  of  0.001  and  the  corrections  due  to  the  wall  constraint  is  thus  negligible  in  the  subsonic  limit. 

At  transonic  speeds  the  effect  of  the  tunnel  constraint  will  depend  on  the  ratio  of  the  transverse  scale 
b/e“ and  the  tunnel  height  and  width  as  discussed. 

For  the  axisymmetric  equivalent  bodies  the  drag  at  zero  lift  is  taken  directly  from  the  ex¬ 
perimental  data  and  no  further  processing  is  needed.  The  cross-sectional  area  of  the  body  B1  is 
identical  to  the  wing-body,  the  blockage  effect  is  very  small  at  the  subsonic  flow  range  and  no  correc¬ 
tion  is  applied.  In  the  transonic  range,  the  outer  flow  field  is  equivalent  to  that  of  the  wing-body, 
thus  the  tunnel  flow  must  be  considered  as  a  complete  system. 

5.0  EXPERIMENTAL  RESULTS  AND  DISCUSSIONS 

5.1  Drag  at  Zero  Lift 

The  experimental  results  are  first  checked  against  the  classical  area  rule.  The  transonic  wave 
drag  is  deduced  from  the  total  drag  by  subtracting  the  viscous  drag  and  the  pressure  drag  which  is  the 
total  drag  before  the  drag-rise.  At  small  and  zero  incidence,  the  drag  in  subcritical  flows  is  dominated 
by  viscous  drag.  For  the  present  experiment  the  drag  at  Mach  number  of  0.8  is  used  as  the  viscous 
drag.  The  transonic  drag-rise  is  shown  in  Figure  7  presented  in  the  similarity  form  as  defined  in 
Equation  (11).  The  drag  parameter  is  normalized  by  the  maximum  cross-sectional  area  SL.  of  the 
model  and  the  transonic  similarity  parameter  K  is  in  terms  of  slenderness  ratio  rb  based  on  Sc.  With 
the  parameters  in  these  forms  the  data  of  the  wing-bodies  and  the  equivalent  bodies  can  be  compared 
directly.  The  drag-rise  of  the  wing-body  WB2  and  the  equivalent  body  B1  are  very  close  through  the 
Mach  number  range  considered.  The  wing-body  WB1  with  a  wing  of  smaller  aspect  ratio  but  a  thicker 
section  gives  persistently  higher  drag  than  the  model  WB2.  This  seems  to  be  contradicting  the  clas¬ 
sical  area  rule,  which  works  better  for  wings  with  small  aspect  ratios*18’22*.  This  discrepancy  is 
believed  to  be  caused  by  stronger  shock  loss  in  the  near  field  and  shock  boundary  layer  interaction 
for  the  WB1  model. 

The  boundary  layer  transition  on  the  wing  is  natural  and  occurs  about  one  inch  (25.4  mm) 
downstream  parallel  to  the  leading  edge  for  both  models.  This  relatively  large  extent  of  laminar 
boundary  layer  interacts  with  the  wing  shock  causing  flow  separation  in  the  wing  tip  region,  can  be 
seen  in  the  flow  visualization  on  the  wing  surface,  Figure  8.  The  separation  appears  to  be  more  severe 
on  the  WB1  model  with  a  thicker  wing.  The  strong  shock  wave  in  the  near  field  of  the  model  WB1 
also  causes  higher  drag.  Drag  due  to  shock  loss  is  proportional  to  [M] 3  where  [  ]  denotes  the 
difference  of  conditions  in  front  and  behind  the  shock  and  in  transonic  flows  this  can  be  shown  to 
be  proportional  to  [U] 3  for  small  disturbances*23*.  Near  the  wings  the  local  drag  due  to  the  shock 
loss  per  unit  spanwise  length  can  be  approximately  related  as  Dj  =*  (1.8)3  D2 ,  where  the  thickness 
ratio  of  the  two  wings  is  1.8.  Integrating  along  the  span  and  noting  that  b,  =  b2  /1.8,  we  have  D!  = 
3.24  D2 .  Thus  in  the  near  field  the  shock  loss  of  the  model  WB1  is  much  higher  than  that  of  WB2, 
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consistent  with  the  observed  drag-rise  as  shown  in  Figure  7.  This  confirms  further  that  the  area  rule 
works  better  with  thin  wings;  for  example,  the  excellent  correlations  obtained  by  Whitcomb  apply 
for  models  with  wing  thickness  of  4  percent  which  generate  relatively  weak  shocks  in  the  near  field. 

The  close  correlation  of  the  drag-rise  of  the  wing-body  and  the  equivalent  body  confirms 
the  equivalence  rule  at  zero  lift.  The  shock  structures  of  the  wing-body  and  the  equivalent  body  are 
shown  in  the  Schlieren  flow  visualizations  in  Figure  9  and  indicate  the  similitude  of  the  outer  flow 
fields.  In  the  near  field  of  the  wing-body,  the  shock  is  formed  near  the  trailing  edge  of  the  wing  and 
is  stronger  for  the  model  WB1  with  a  thicker  wing  than  the  model  WB2.  For  the  equivalent  body,  the 
shock  is  located  at  the  rear  part  of  the  bulge  and  is  more  oblique  near  the  body  surface.  The  boundary 
layer  does  not  appear  to  separate  at  the  rear  end  of  the  bulge,  nor  strong  shock-boundary  layer  inter¬ 
action  is  indicated.  The  shocks  on  all  models  change  rapidly  to  axisymmetric  normal  shocks  within  a 
transverse  distance  of  two  body  diameters,  forming  closely  similar  outer  flow  fields. 

The  outer  flow  similarity  can  also  be  examined  from  the  pressure  distributions  along  the 
upper  and  lower  walls.  Figure  10  shows  the  variation  of  the  pressure  coefficient  along  the  walls  for  all 
three  models  at  Mach  number  of  0.98.  Similarity  of  the  outer  flow  fields  of  all  three  models  are 
observed.  The  outer  flow  pressure  variations  reflect  directly  the  source  distribution  along  the  body 
axis  and  only  weakly  retain  the  rapid  variation  of  the  flow  quantities  near  the  body.  At  other  Mach 
numbers,  the  outer  flow  similitude  has  also  been  observed  from  the  Schlieren  flow  visualizations 
and  the  pressure  distributions  near  the  walls.  The  trend  is  closely  similar  to  the  cases  just  discussed. 


5.2  Lift 


The  lift  vs  angles  of  attack  curves  for  both  models  have  already  been  presented  in  Figures  5 
and  6.  The  linear  portion  ofthe  lift  curve  before  onset  of  flow  separation  can  be  represented  by  the 
lift  curve  slope  at  zero  lift,  CL  .  In  Figure  11,  the  lift  slopes  are  correlated  in  the  form  of  the  clas- 
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sical  transonic  similarity  law,  which  applies  to  the  present  models  with  geometric  similar  wing  plan 
forms  and  sections.  At  subsonic  speeds  both  models  follow  a  similar  trend  with  increasing  lift  slope 
as  Mach  number  increases  but  fail  to  reach  the  theoretical  value  of  7r/2  at  Mach  unity.  At  supersonic 
speeds,  the  correlations  are  less  coherent  and  the  lift  slopes  decrease^ as  Mach  number  increases.  At  the 
low  Mach  number  end,  the  corresponding  incompressible  values  of  C,  and  A  obtained  from  applying 

the  linear  compressibility  transformation  follow  the  existing  data  for  delta  wing  and  the  lifting  surface 
theory  closely*24*.  At  Mach  number  unity,  the  values  follow  the  transonic  correlation  in  the  form 
Ci  (t/c)*/3  vs  A(t/c)'/3  *25*.  The  model  WB1  with  smaller  aspect  ratio  is  within  the  linear  range 

CK 

(slender  body  theory),  while  the  model  WB2  with  larger  span  is  beyond  the  linear  range. 

As  the  flow  Mach  number  approaching  unity,  the  lift  slope  Cl  depends  strongly  on  the 

shock  movements  at  the  upper  and  the  lower  surface  of  the  wing*2*  ,27) .  The  large  dip  at  M^  =  0.96 
for  the  model  WB1  is  caused  by  the  movement  of  the  lower  surface  shock  towards  the  trailing  edge 
while  the  upper  surface  shock  remains  nearly  stationary  at  small  positive  incidence.  This  increases 
the  total  suction  at  the  lower  surface  and  hence  reduces  the  lift.  At  higher  Mach  number,  both  shocks 
move  to  the  trailing  edge  and  the  lift  recovers.  For  the  thin  wing,  the  difference  of  the  shock  move¬ 
ment  is  less  severe  and  loss  of  lift  does  not  occur. 

When  the  angle  of  attack  increases  to  about  the  value  of  the  wing  thickness,  the  flow 
separates  from  the  leading  edge  and  rolls  up  into  a  vortex.  This  can  be  seen  clearly  from  the  surface 
streamline  patterns  as  shown  in  the  flow  visualization  pictures  Figure  12.  For  the  model  WB1,  at 
5  degree  incidence,  the  leading  edge  vortex  is  still  weak  and  lies  close  to  the  leading  edge.  At  10  degree 
incidence  a  strong  conical  vortex  is  located  along  a  generator  of  the  delta  wing  well  in-board  from  the 
leading  edge.  A  similar  flow  pattern  can  also  be  observed  over  the  wing  of  the  model  WB2  at  inci¬ 
dences  about  one  half  of  that  of  WB1.  The  lift  induced  by  the  vortex,  however,  is  less  significant  than 
that  on  the  highly  slender  wings*28 ,29*.  The  lift  curves  do  not  appear  to  have  appreciable  nonlinear 
lift  contributions  even  at  subsonic  Mach  numbers  (see  Figs.  5  and  6). 


The  equivalence  rule  outlined  in  Section  2  are  based  on  the  solutions  of  the  attached  flow*5! . 
Flows  with  leading  edge  separation  has  been  treated  by  Barnwell*6  7  *  and  the  results  are  closely 
similar  to  those  given  in  Reference  5.  Thus  the  equivalence  rule  in  the  present  form  will  apply  to  the 
complete  range  of  CL  investigated  in  the  experiment. 

5.3  Induced  Drag 


The  drag  measured  in  the  experiment  consists  of  contributions  from  the  skin  friction, 
pressure  forces  including  those  due  to  shock  waves  and  drag  induced  by  lift.  At  subcritical  conditions, 
there  are  no  shock  waves  and  the  pressure  drag  is  much  smaller  that  the  viscous  drag.  Thus  before  the 
transonic  drag  rise  the  total  drag  contains  the  viscous  drag  and  the  lift  induced  'frag  only.  At  zero  lift 
the  induced  drag  is  zero  and  the  measured  drag  can  be  considered  as  mainly  skin  friction  drag.  The 
boundary  layer  flows,  however,  depend  on  the  incidence  and  as  the  Mach  number  increases,  shock 
boundary  layer  interaction  may  also  occur,  thus  the  viscous  drag  will  depend  on  both  the  flow  Mach 
number  and  the  angle  of  attack  of  the  model.  However,  the  effect  of  incidence  on  the  viscous  drag, 
before  the  onset  of  separation  on  the  wing  surface  is  known  to  be  small.  Thus  for  the  present  study 
the  skin  friction  drag  is  assumed  to  be  invariant  and  have  the  value  equivalent  to  the  drag  measured 
at  Mach  number  of  0.8  and  zero  angle  of  attack. 

The  drag  induced  by  lift  is  well  established  in  the  wing  theory*25'30'.  For  a  delta  wing 
with  small  aspect  ratio,  the  slender  wing  theory  gives  the  induced  drag 


Cn  = 


C£ 

7T  A 


(27) 


which  is  the  optimized  value  corresponding  to  an  elliptical  loading.  For  a  genera]  planform, 


*1 


Gp.  =  k  — 
7T  A 


(28) 


where  k  is  the  induced  drag  factor  which  depends  on  the  planform  shape.  The  induced  drag  can  also 
be  considered  as  a  component  of  the  normal  force  in  the  free-stream  direction  subtracting  the  suction 
force  that  may  exist  at  the  leading  edge  of  the  wing.  If  there  is  no  suction  force,  the  induced  drag  is 
simply 


'D; 


=  C,  a 


(29) 


which  is  twice  as  large  as  that  when  a  full  suction  force  is  present,  Equation  (27).  For  a  general  plan- 
form,  the  induced  drag  falls  in  between  these  two  extremes  and  is  a  function  of  the  lift  and  the  lift 
curve  slope*25).  It  is  written  in  a  general  form  as 


Cn,  =  k 


(30) 


where  k  is  the  induced  drag  factor. 


The  drag  polars  for  both  models  at  Mach  number  0.8  are  replotted  in  terms  in  Figure  13 
for  both  models.  The  limits  for  a  full  suction  force  and  no  suction  are  also  shown.  The  experimental 
data  fall  in  between  these  limits.  At  low  lift,  the  leading  edge  suction  reaches  about  a  half  of  the  full 
value  and  decreases  steadily  as  the  lift  increases.  The  correlation  also  shows  that  CD  is  not  a  linear 
function  of  C  2 .  Thus  the  general  relation  Equation  (30)  should  be  used  to  correlate  the  induced  drag 
with  lift. 
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As  the  Mach  number  increases  the  drag  due  to  formation  of  shock  wave  appears.  For  the 
wing-body  configurations  there  are  two  shocks  formed  over  the  model,  one  at  the  end  of  the  nose 
section  and  the  other  at  the  rear  of  the  wing.  From  the  experimental  data  of  the  axisymmetric  bodies, 
it  is  found  that  the  variation  of  the  nose  wave  drag  with  incidence  is  very  small.  Thus  the  increase  of 
drag  with  incidence  is  attributed  mainly  to  the  wave  drag  and  the  induced  drag  due  to  lift  of  the  wing. 
The  drag  data  with  the  zero  lift  value  CD  subtracted  are  plotted  in  terms  CdCLq  against  C2  in 

Figures  14  and  15  for  the  two  wing  body  models  respectively.  Most  of  the  data  are  successfully  cor¬ 
related  by  a  straight  line.  At  the  high  lift  end,  the  data  deviate  from  the  straight  line  correlation  as  the 
values  of  CL  decrease  rapidly  at  high  incidences  with  flow  separation  occurring  on  the  upper  surface 
of  the  wing. 

5.4  Lift-Depending  Wave  Drag 

The  slope  of  the  correlation  curve  discussed  in  the  preceding  section  gives  the  overall  lift- 
depending  drag  factor  which  consists  of  the  induced  drag  due  to  lift  and  the  lift-depending  wave 
drag<31).  The  drag  induced  by  lift  is  independent  of  Mach  number^24*.  Thus  the  value  of  k  at  Mach 
number  of  0.8  before  the  drag-rise  can  be  taken  as  the  induced  drag  factor  for  all  cases.  The  lift- 
depending  wave  drag  can  now  be  deduced  as 


(k  ^0.8 ) 


C2 


The  values  of  k  as  a  function  of  Mach  number  are  shown  in  Figure  16. 


(31) 


The  lift-depending  wave  drag  is  now  presented  in  terms  of  the  parameters  of  the  equivalence 
rule.  The  conditions  of  equivalence  require  the  doublet  and  the  effective  source  to  be  identical  at  the 
same  transonic  similarity  parameter  K.  From  Equation  (24),  the  parameter  controlling  the  lift  is 
e'ai  ff  which  can  be  written  in  terms  of  the  experimental  variables  as 


a„ 


(7  +  1)*/2  co  1 


8  (X'r),/2 


(M„ACl) 


(32) 


For  the  present  experiment,  (X  t),  cq  /C  are  identical  for  both  models,  thus  the  lift  controlling  param¬ 
eter  can  be  expressed  as  MwTCL .  The  drag  coefficient  is  now  referenced  to  the  maximal  cross-sectional 
area  of  the  model.  Thejesults  in  terms  of  these  parameters  are  shown  in  Figure  17.  The  abscissa  is 
given  in  two  scales,  M^X  CL  as  in  Equation  (32)  and  e‘av/T'  =  or*  |  Cn  e'|"l/2  as  >n  the  original  form  of 
the  parameter,  Equations  (3)  and  (24).  The  results  are  presented  from  Mach  number  of  0.9  to  1.1. 
Below  0.9  the  wave  drag  is  negligibly  small.  The  values  of  the  transonic  similarity  parameter  K  is  given 
beside  the  Mach  number  for  each  case.  For  most  cases,  the  correlation  of  the  wave  drag  with  the  lift 
parameter  is  found  to  be  similar  for  both  models.  The  correlation  demonstrates  that  for  the  models 
with  the  same  cross-sectional  area  but  different  aspect  ratios  and  thicknesses  of  the  wings,  if  the  lift 
and  the  effective  cross-sectional  area  due  to  lift  are  the  same  at  a  given  transonic  similarity  parameter, 
then  the  wave  drag  due  to  lift  is  the  same.  The  equivalence  rule  with  lift  is  thus  established  experi¬ 
mentally.  For  the  case  of  Mach  number  0.96,  the  model  WB1  has  a  much  larger  drag-rise  due  to  lift 
than  other  cases.  This  is  caused  by  the  low  value  of  the  lift  curve  slope  at  this  particular  Mach  number 
as  discussed  in  Section  5.2.  Thus  in  order  to  achieve  the  same  lift  condition  as  defined  by  the  lift 
parameter,  the  model  has  to  be  at  much  higher  incidence  than  in  the  other  cases  and  hence  larger  drag 
is  obtained. 


For  the  lifting  condition  with  the  lift  parameter  e  ac  It  of  0.2,  the  equivalent  body,  B2, 
with  effective  area  due  to  lift  added  to  the  basic  body,  Bl,  was  tested  in  the  same  Mach  number  range 
(see  Section  3.5).  The  drag-rise  of  the  equivalent  body  and  the  wing-bodies  at  the  design  lift  condition 
are  shown  in  Figure  18  presented  in  terms  of  the  similarity  parameters.  The  difference  in  zero  lift 
drag  of  the  wing-bodies  as  shown  in  Figure  7  is  again  reflected  in  the  higher  drag-rise  for  the  WB1 
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model.  The  drag-rise  of  the  equivalent  body  B2  is  slightly  lower  than  that  of  the  wing-body  WB2.  This 
discrepancy  is  due  to  the  boundary  layer  separation  at  the  rear  of  the  hump  of  the  equivalent  body 
giving  an  effective  fairing  to  the  body  contour,  much  less  abrupt  than  the  physical  one.  This  effective 
broadening  of  the  hump  generates  a  weaker  shock  than  that  which  may  have  formed  if  there  had  been 
no  separation  and  hence  a  lower  drag-rise.  The  boundary  layer  separation  at  the  rear  of  the  hump  and 
the  reattachment  downstream  at  the  aft  body  can  be  seen  clearly  from  the  surface  oil  flow  visualiza¬ 
tions  in  Figure  19.  The  rounding  off  of  the  wing  tips  to  produce  a  more  gradual  change  of  the  body 
contour  as  discussed  in  Section  3.3  does  not  seem  to  be  adequate  in  preventing  flow  separation. 

The  Schlieren  flow  visualizations  for  both  wing-body  and  the  equivalent  body  models  at  the 
design  lift  condition  are  shown  in  Figure  20.  Since  the  Schlieren  pictures  are  taken  at  a  fixed  time 
interval  during  continuous  pitching  sequence  of  the  test,  the  lift  parameter  cannot  be  matched  exactly 
for  both  models  and  the  condition  shown  in  the  Figures  are  the  ones  closest  to  the  design  lift  condi¬ 
tion.  The  shock  patterns  over  the  wings  of  both  wing-body  models  are  closely  similar.  It  should  be 
noted  that  for  the  same  lift  parameter  the  angle  of  attack  is  not  the  same  for  the  models.  The  shock 
pattern  over  the  rear  part  of  the  body  B2  is  more  complicated  than  those  of  the  wing-bodies.  It 
indicates  clearly  a  strong  interaction  of  the  shock  with  the  boundary  layer  separating  at  the  rear  of 
the  hump  of  the  body. 

The  similarity  of  the  far  field  can  be  further  examined  by  the  static  pressure  variations  near 
the  tunnel  walls  as  in  the  zero  lift  cases.  The  static  pressure  distirbutions  at  Mach  numbers  of  0.94 
are  shown  in  Figure  21.  The  pressure  coefficient  is  now  written  in  the  similarity  form  as  defined  in 
Equation  (10),  so  that  the  data  of  the  equivalent  body  B2  which  has  a  slightly  larger  maximal  cross- 
section  can  be  compared  in  the  same  scale.  Because  the  pressure  scanning  was  done  at  predetermined 
incidences  as  described  in  Section  4.1,  the  lift  parameters  of  two  models  cannot  be  matched  exactly. 
The  cases  shown  in  the  Figures  are  the  closest  matches  that  can  be  obtained,  although  the  lift  param¬ 
eter  has  a  larger  value  than  the  design  condition  for  the  equivalent  body  B2.  Taking  into  consideration 
the  data  scattering  and  the  pressure  signal  being  very  weak,  the  data  for  both  models  follow  the  same 
trend  closely.  At  lifting  conditions  the  upper  and  the  lower  wall  pressure  distributions  are  no  longer 
the  same.  At  the  lower  wall  a  strong  compression  directly  underneath  the  models  can  be  observed. 

The  pressure  distributions  for  the  equivalent  body  B2  are  also  plotted  in  the  Figure.  Because  of  the 
difference  in  the  lift  parameter,  the  data  cannot  compare  quantitatively  with  the  wing-body  results. 
Qualitatively,  the  data  at  the  upper  wall  follows  the  trend  of  the  wing-body  data  closely,  since  the 
lifting  wing  and  the  humped  body  both  have  the  same  flow  pattern  of  an  expansion  followed  by  a 
recompression  shock.  Having  a  smaller  lift  parameter,  the  magnitude  of  the  B2  data  is  consistently 
lower  than  those  of  the  wing-body.  The  B2  data  at  the  lower  wall  is  closely  similar  to  those  of  the 
upper  wall  because  of  axial  symmetry  and  do  not  have  a  region  of  compression  as  generated  by  a 
lifting  wing. 


6.0  CONCLUSIONS 

The  transonic  equivalence  rule  with  lift  has  been  studied  experimentally.  The  experiment 
has  established  the  similitude  of  the  outer  flow  fields  of  two  wing-bodies  under  the  conditions  defined 
by  the  equivalence  rule(5).  The  two  models  have  the  same  cross-sectional  area  distribution  but  with 
wings  of  different  aspect  ratios  and  thicknesses.  Two  equivalent  bodies  of  revolution  have  also  been 
tested,  one  corresponding  to  the  zero  lift  condition  and  the  other  with  lift.  The  experimental  results 
are  summarized  as  follows. 

The  equivalence  rule  at  zero  lift  works  well  for  the  wing-body  model  with  a  thin  wing  with 
thickness-chord  ratio  0.044.  For  the  wing  with  thickness  of  8  percent,  the  rear  field  drag  caused  by 
the  strong  shock  over  the  wing  gives  a  higher  drag-rise  than  that  of  the  equivalent  body. 

The  correlations  of  the  lift-depending  wave  drag  and  the  lift  parameter  e  at  /r  of  two  models 
are  closely  similar  for  the  tested  Mach  numbers.  Similitude  of  the  outer  flow  field  under  the  same  lift 
condition  defined  by  the  lift  parameter  is  established. 
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The  drag-rise  of  the  equivalent  body  with  lift  correlates  reasonably  well  with  that  of  the 
wing-body  at  the  design  lifting  condition,  though  the  drag  parameter  has  a  slightly  lower  value 
throughout  the  Mach  number  range.  This  is  caused  by  the  flow  separation  at  the  rear  of  the  hump  of 
the  body  forming  an  effective  contour  which  reduces  the  local  shock  strength. 

The  experiment  therefore  verifies  the  conditions  defined  by  the  equivalence  rule  with  lift. 

It  establishes  the  similitude  of  the  outer  flow  field  at  lifting  conditions.  As  in  the  classical  rule  of  zero 
lift,  the  equivalence  rule  with  lift  works  well  for  configurations  with  a  thin  wing  and  a  gradual  varia¬ 
tion  of  the  equivalent  cross-sectional  area  which  includes  the  geometric  one  and  the  effective  area 
from  the  axial  lift  distribution. 
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TABLE  1 


body  length 
body  radius 
root  chord 
semi-span 
wing  thickness 
aspect  ratio 


WING  SECTIONAL  PROFILE 


x/c 

t/c 

0 

0.00000 

0.05 

0.02700 

0.10 

0.03573 

0.20 

0.04502 

0.30 

0.04914 

0.40 

0.04989 

0.50 

0.04776 

0.60 

0.04284 

0.70 

0.03519 

0.80 

0.02498 

0.90 

0.01234 

1.00 

0.00000 

TABLE  2 

GEOMETRIES  OF  THE  MODELS 


MODEL  WB1 

Cj  =  17.5  in.  (444.5  mm) 
rj  =  0.875  in.  (22.23  mm) 

cq  =  6  in.  (152.4  mm) 

b,  =  4.0437  in.  (102.71  mm) 

tj  /c  =  0.8 

A  =  2.60 

X,  =  0.231 

r,  =  0.0516 


MODEL  WB2 

C2  =  17.5  in.  (444.5  mm) 
r2  =  0.875  in.  (22.23  mm) 

cp  =  6  in.  (152.4  mm) 

b2  =  7.2856  in.  (185.05  mm 

t2/c  =  0.4444 

A  =  4.69 

X2  =  0.416 
t2  =  0.02865 
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TABLE  3 


EFFECTIVE  AND  GEOMETRIC  CROSS-SECTIONAL  AREAS  AND 
RADII  OF  THE  EQUIVALENT  BODIES  OF  REVOLUTION 


X'Xo/Co 

Sc 

SLift 

Se 

rc  (in) 

re  (in) 

0.00 

0.6632 

0.0000 

0.6632 

0.8750 

0.8750 

0.50 

0.6632 

0.0012 

0.6644 

0.8750 

0.8758 

0.10 

0.6632 

0.0043 

0.6675 

0.8750 

0.8778 

0.15 

0.6728 

0.0091 

0.6819 

0.8836 

0.8872 

0.20 

0.7020 

0.0155 

0.7175 

0.9026 

0.9101 

0.25 

0.7369 

0.0234 

0.7603 

0.9254 

0.9369 

0.30 

0.7764 

0.0328 

0.8092 

0.9492 

0.9665 

0.35 

0.8161 

0.0436 

0.8597 

0.9732 

0.9962 

0.40 

0.8564 

0.0558 

0.9122 

0.9968 

1.0262 

0.45 

0.8937 

0.0693 

0.9630 

1.0184 

1.0544 

0.50 

0.9287 

0.0841 

1.0128 

1.0380 

1.0813 

0.55 

0.9578 

0.1001 

1.0579 

1.0542 

1.1051 

0.60 

0.9813 

0.1174 

1.0987 

1.0672 

1.1262 

0.65 

0.9956 

0.1359 

1.1315 

1.0748 

1.1429 

0.70 

1.0000 

0.1556 

1.1556 

1.0774 

1.1550 

0.75 

0.9917 

0.1765 

1.1682 

1.0728 

1.1613 

0.80 

0.9685 

0.1985 

1.1670 

1.0600 

1.1607 

0.85 

0.9262 

0.2217 

1.1479 

1.0368 

1.1512 

0.90 

0.8621 

0.2460 

1.1081 

1.0002 

1.1310 

0.95 

0.7704 

0.1570 

0.9274 

0.9456 

1.0347 

1.00 

0.6632 

0.0471 

0.7283 

0.8750 

0.9055 

k, 
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TABLE  4 

AERODYNAMIC  COEFFICIENTS  OF  WING-BODY  MODEL  WB1 


MACH  NO 
R E/FT ( M 

A I  PH 8 B 

-4.6104 
-4.  o  1R6 
-3 .4654 
-2.9  411 
-2. 421  C 
-1.9565 
-1  .  3268 
-0 . 8 1 5  R 
-0.2  06  6 
u  .25?” 

0.  769i 

1 .8  61  R 

">  1  QT  •} 

2*.  9  331 

3  .  4  6''  6 

4  .  C  1  4  9 
4.  5427 
5.0690 

5  .  6  90  R 

6.  V31  7 
6 . 6  64  1 
7.24^9 

7. 

8.2931 
8. 822'' 

9 . 353  1 
o.  9  77  4 
10. 4122 
12.  5  50  0 


)  =7.  Q6 

CL  R 

-5. 2 7 56 
. 2  r  R  5 
-7 .1  R3  5 
-3. 1 Rf 7 
28 

— ' '. 1  ’  4  0 
-”  .C  546 
-')  .957  5 

-c.o  in 

O.i  1  ?q 

5  .  5  ?  1  1 

'.'625 
2 . : 94 R 
r . 1 19  3 
0 . 1  4  6 1 
''.1699 

o  .  2  ”  R  6 
0.2276 
7.?6U 7 
”  .  2  7  R  9 

''.333'’ 

C  .  32qo 
5  .  7  7  2  5 
T.  ?d(iq 
0 . 4  ■*  6  B 

r .  4  072 
0.4618 
5.4767 
0.5061 
n  .  6  1 4  6 


0.5  27 ? 
~  .  c  29  1 
5  .  r  1  -  9 
:  .  r  1 5  8 
2  .  C  1 4  2 
”  .  r  1 2  o 
;*■ .  r  1 1  9 

5 . 5 1 1  6 

2.011  1 
'',0111 
■’  .  C  1  i  5 
”.712? 

5.0138 
”  .  ',1C  6 
9  .  r  1 C,  8 
1.OI04 
7  .  C  2  1  0 
7.  ”24  9 
0  .  C  29  4 
C  .  7  37  3 
9 . 0374 
0.044-’ 

5 .  c  5 1 2 

9. 0  57  1 

3.0642 
”  .  07i 4 
7 .070  2 
7.0990 
7.0921 


MACH  NO. -7, 714 
RE/*? (Ml  =9.31 
AI.PHAB  PL  9 


-4. 7785 
-4.  i'-<--7 
— 3 .  5  84  7 
-3 . 0426 
-2. 521 7 

—  7.9  62  4 
-1.41.7 
-0.8596 
-0.3  46  4 
0.1917 

a  7  £j  ->  -i 

1 1  292  4 

1.  9529 

2 . 3  597 

2.  9m 
3.4264 
7.  971  P 

4.  524  7 

5.  ['59  4 
5 . 6022 
6.1417 
6.6825 
7.2467 
7. 766,9 
8. 2946 
8.8211 
0. 3524 

9.870? 
1r  .  4196 
10. 561  6 


-7.267  3 
-''.2  22  4 
-0.1927 
-3.17”  /> 
•>  r  6 

-o:cei2 

-3.0765 
-r> .  -  a  5 

i .  o :  6 1 
A . '  307 
0 .  o  7  6  rt 
5 . 5  P  9  7 
A.72i 8 
7. 164  1 
0.173  9 
?  .  7  ”  o  7 
0 . 2  4  7  3 
7.2609 
O  ,  3  'T  7 

9 , 9  ?'■>  5 

9.3560 

7.3944 

7.4-11 

9.4233 

7 .4  UOq 

7.4  66  2 

1.4  R6  8 

9.5  156 

9-.61Q9 


”  .  r op  1 

7 .0  22'' 
”  .  n  1  q  2 
A  .  C  1 7  1 
7.3149 

1/'V) 

7.9127 

7.0174 

7.5-774 

7.511  7 
7.0116 

0 . 7  1 3  1 
3 .5 1 37 

7.^14  7 

7.0160 
A  .  0  1  7  6 
0. 0  2:  4 

7  .  ”  7  3  Q 

7.7278 

r.  q  - 1  -> 

7*C  36  6 
7,3479 
7  .  c  57  7 
9 .756  5 
■'.  0  622 
0.57-6 
7.0775 
r  .''86  9 
0  .  r  Q  6  2 
■:  .003  7 


MACH  NO. 
RE/  FT  (M) 

ALPHAB 

-4.77P1 
-4. 1 226 
-3.  5  01  7 
-3 . 0  53  6 
-2. 5  31  9 
-1  .  99-36 
-1 .4170 

-0. 90-50 
-5.  3-2  5 
0.1 847 
7,7111 
1.2487 
1  .  7982 
2.3299 

2.  8631 
3.4377 

3.  9487 

4. U9GC 
5 . 017? 
5.5771 

6.H61 
6.6363 

7.  1523 
7.6930 

8. 2477 
8.7711 
9 , 3-73 
9.8256 
17. 3768 
1C. 5367 


•  0 . 8  4  o 

=9.02’ 


-7.2652 

-0.227  0 
-7, 1 Rc  9 
-0.1 673 
-1.1  744 

-".7765 

-0 . r f  a  1 
-0.r  57  3 
-''.--iqi 

C  .  C  C  8  8 
7.7373 

A.3644 
C  . ”  04 r 
9.1796 
O.1404 
7.1941 

0.2rQ2 
0.2 356 
'.26  17 
3 .  ?072 

7.3  26.3 
■5 .  7  4  o  2 
0.366  1 
3.7  97  3 
0.4149 
3.4418 
7.4586 

3.4  85  7 
O . 6  QO  R 

0.5161 


7.0251 
o. ”2i 4 
■9.0104 
5- .  c  1 6  2 
.0149 

7.5173 

r  ^  i  -  q 

.‘OUR 
”, r  1  -  4 
5.7114 
7  .  C  11  6 

,0  p  .  -V  , . 

0*5129 
".513° 
5.0166 
7 . 0  1  7  ” 
7.0198 
7.0296 
3  .  C  26  1 
5 .73-  3 
'9.7  356 

C  .  7  a”.  7 
7 .7  457 
0.0  52  4 
3  .  C  5o  4 
7.75-2 
7.0747 
0.5933 
7 .093  5 
0.096  3 


MACH  NO. 

3E/FT ( M ) 

AL  ?HA  P 

-4.713  4 
-4. 1 37- 
-3.  5  63  6 
-3. ”473 
-2.5  39? 

-1 . 941 5 
-1 . 4275 
-3.3  999 
-C.3330 
0.2017 
9. 7468 
1.2757 
1.9146 
2.3692 

2.  9935 

3 . 4  77  7 

3. a047 

4 . 5  4  7  4 

5. ” 957 
5.5939 
6. 1738 
6.7241 
7.  2  51  a 
7.7521 
8. 3058 
3.8466 
0. 3585 
9.8879 

r  .  U51  8 

10. 5  91  3 


=  7.89  0 

=  9.  04 
CL  B 

-9 .2644 
-0.2331 

-9.202  7 

-i.iqri 

-9.1467 

-3.1043 

-3.3772 

-5.75m 

-3.5181 

7.7152 

7.9385 

9.7cag 
7  .  ■”  9  9  4 
7.1  ;,"R 

0.1621 

■3 . 7  6  5 

0.2224 

0.2611 

7.2805 

0.2988 

7.3448 

0.37-11 

7. 3082 

3 . 4  C  6  Q 

7.4  208 
3.4547 
0.4  15  1 

0.4954 

7.5257 

0.53C7 


3.0  qi 
5.0238 
O.cif 6 
”  .  r  3  7  a 
3.0167 
0.0136 
0.0126 
C  .0127 

3.5116 

0.0116 

0.0119 

7.0723 

0.5132 

7.7149 

0.0157 

0.5196 

0.7  221 
5 . 7*26  1 
7.033  5 
0  .  C  3  7  9 
0 .049  7 
7.546  9 
0.5532 
C  .  0590 
0.0655 
0.7717 
7/ .  C  8 1  4 
0.0897 
0.105  5 
0. 1028 


1 
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TABLE  4 

AERODYNAMIC  COEFFICIENTS  OF  WING-BODY  MODEL  WBI  (Cont’d) 


MACH  no.' 
7E/FT(M)  ' 
A  T.  P  HA  H 

-4.  6  644 
-4.  0  70  2 
-3.  5  358 
-2.9895 
-2.501'’ 

-1 . 9935 
-1. 371 r 
-0.3293 
-0.23  98 
C.2  522 
" . 9^30 
1 .3490 
1 . 3336 
2.4330 

2.  9  629 
3.5355 
4.0616 
4 .6  307 
5.  1  44  4 
5.6954 
5.  2207 
6.7544 
7.  302  6 
7.8276 
0. 3623 

8.8  98  7 

q. 4238 
9. 9583 
12. 4  94  9 

10.6532 


=  7.918 
=  3.  03 
CL  B 

-7.2665 
-ft . 2  25  3 
-9.2007 
-0.1794 
-9.1 322 

-o.  i  -->4 

-3.0851 
-0.  '.'4C  7 
-2. r 120 
9.0113 
O.-'u0!? 

3.0779 

n.i120 

■1,1410 

7. 1675 

9.2057 

0.2348 

0.2  58  5 

9 . 2  °9  2 

3.3249 

7.3468 

0.3  682 

0.3  93  7 
3.4178 
0.439  9 

7.4616 
9.478H 
9.593 1 
7.5272 
0.5361 


M=  30 


CD 


0.0  30  7 

0.0257 
3.C22  1 
3.0192 

7.9171 

0.0145 
3 . 0  1 1 1 
7.3126 
9.0121 
3.7122 
0.0124 
0.0132 
<3.0143 
0.7162 

7. 7  ifi  3 
0.7715 

0.0249 

0.0288 

9.7  332 

0.9389 
0.0  43  5 

0 . 049  7 

0.0  55  3 
0. 0623 
0.0694 
0  .077  C 
■7.9843 
1.0933 
0.  102  4 
7.106? 


JJACH  NO. =0.939 

OE/FTH)  =8.  CO 


alpha? 

-4. 7017 
-4.  1  59  1 
-3 .5819 
-3.0432 
-2.5494 
-1.976? 
-1.4374 
-0.9122 
-0 .3  339 
0.2255 
0. 7376 
1.2738 
1.8323 
2.3699 
2.9193 
3.4676 
3.9979 
4.5379 
5. 0383 

5.  6260 
6. 1838 
6.721  5 
7. 2555 
7.7989 
8. 3312 

8,8  73  8 
9.  3957 
9.  9190 
17. 4699 
10.6312 


CL  R 

-0.2546 
-0.2353 
-0.2042 
-0.1789 
-0.V557 
-0.1261 
-0.0 8? 4 
-0.0  618 
-0.0 166 
0.0319 
0.CSC5 
9.0834 
7 . 1  164 
9.1441 
0. 1 709 
7.1943 
7.2258 
0.2530 
0.2869 
9.3195 
0.3470 
9.3845 
0.4076 
0.4240 
0.4455 
0.4772 
7.4877 
9.5166 
0.5382 
0.5446 


1=30 

CD 

0.0340 
7 . 0  3C  1 
0.0258 
C. 0228 

0,0  Of'  5 

0.C18  1 
C.016  1 
0.0152 
0.0145 
0.0145 
0  ,  C  1  4  9 
0.0159 
0.0173 
0.0192 
7. 0  21  7 
0.024 5 
0  .  C  28  2 
0.0320 
0 .0  368 
0.0421 
0.7477 
0.0546 
0.0599 
0.0686 
0.0734 
9.0829 
0.0890 
0.C98  3 
0.  1076 
9. 1  in  4 


MACH  NO.  =  \o<;’ 
85/FT  f 1)  =8. ft 8 ' 
ALPHAS  Cl  B 


8=3' 


CD 


-4.6517 
-4.0  59  2 
-3.5491 
-2.9910 
-2 . 4  76  0 
-1.923? 
-1.3  999 
-0.8339 
-7.3112 
0.2  505 
0. 7792 
1.3491 
1.8757 
2.4079 
2.9412 
3.4913 
4.0188 
4.5804 
5.  1182 
5.673C 
6. 2C33 
6.7596 
7.  3029 
7.8601 
8.3655 
8. 9419 
9.4788 
10.0127 
10. 5595 
10.7191 


-0.2409 
-7.2062 
-ft. 1 01  9 
-9.1 554 
-0.1 316 
-0.17*4 
-0.9839 
-0 . 0  3  9  8 
-7.021 4 
0 .  ft  1 59 
0.C48  1 
9.7883 

ft.  •)  f’p 3 
0.1  30  4 
0.14Q5 

9. 1763 
0.2017 
0.2377 
9.2673 

9.3* 65 

0.3325 
C.3  7C5 
7.408  1 
7.4  36  8 
9.4530 
0.4943 
0.5276 
C .  5442 
0.566  1 
0.5791 


9.0388 

A.0344 

0.0311 

0.02P4 

0.02*2 

C.O747 

9.0233 
C  .  0  22  3 
0 .021 0 

,'.022i 
0. 0223 
0.0234 
0.0242 
7 . 025  9 
0 .0280 
0.030  4 
0 . 0  33  1 
0. 0372 
0.041  5 
0.0472 
ft. 0521 
0.0587 
0.0659 
0,072° 

0 . 079  1 
0.0893 
9.098C 
ft  .  1  06  7 
0.1157 
C  .  1 1 9  9 


JJ AC H  NO, =7. Q3  1 

R5/FTM)  =8.02 
ALPHAB  CLB 


-'1=37 


CD 


-4.7614 
-4. 2092 
-3.6193 
-3. 084  5 
-2 . 5690 
-1. 9929 
-1.4573 
-0.9136 
-0.3567 
0.1773 
0.7598 
1.2749 
1.3253 
2. 3751 
2.910° 
3.4  546 
4.0059 
4.  5?53 
5.1299 
5.6646 
6. 192C 

6.  7509 
7.2777 

7.  8597 
8.3630 
8.9338 
9.4575 
9.  °  990 

10.5519 

10. 6699 


-0 . 7  3Q4 

-9. 2*45 
-9.2046 
-9.191 1 
-0.1552 
-9.1183 
-9.0909 
-9.C544 
-9.9216 
9.9080 
9.0560 
0.0740 
9.1017 
9.1 449 
9.1734 
0.2054 
0.2436 
9.2611 
ft. 3013 
9.344° 
9.3677 
0 . 3  °9  1 
9.4293 
9.4611 
9.4874 
9.5176 
0.5358 
9.5728 
9.6089 
0.6116 


C  .051  5 
9.0 470 
9.0422 
9 .  ?47  2 
9.0378 
9.0355 
9.0342 
C.C332 
0.9326 
0.0324 
0.033  1 
9.C34C 
•9.C35C 
ft. 0370 
9.03°1 
•9.0  41  5 
0.045C 
9.0477 
9.0528 
0.0579 
0.0626 
9.0688 
9,0749 
0.0  82  7 
9.0897 
0.0986 
9.1060 

Ol 1 27  i 
0  .  1 29  1 
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TABLE  4 

AERODYNAMIC  COEFFICIENTS  OF  WING-BODY  MODEL  WB1  (Cont’d) 


MACH  NO. 

=  *1. 99° 

M  =  ?9 

MACH  NO. 

=  1.  rS  1 

<1 

O 

RE/  FT  { M) 

=  7.9  6 

PE/ FT  (Ml 

=p.ci 

ALPHAB 

CL  B 

CD 

ALPHAB 

CL  B 

CD 

-4. 8  32  3 

-3.3070 

9.054  8 

-4. 72r,0 

-0.2755 

9.0600 

-4. 2140 

-3. 2  6  80 

0.0499 

-4.1507 

-U.2494 

3.C560 

-3.6644 

-0.2  23  7 

*  .0454 

-3.6275 

-3.220  1 

0  .  C  52  6 

-3.1 30C 

-9. 1Q!i1 

0.042  4 

-3.0464 

-3.1 88° 

9 . 9  4  0  r 

-2. 6022 

-2.1566 

0.C399 

-2. 5324 

-O.1593 

9  *  fair 

-2.0332 

-0.1251 

0.0378 

-1.9589 

-0.1 381 

C  .044  6 

-1.4935 

-9.1013 

3.0369 

-1.4072 

-0.3739 

9.G491 

-0 . 9692 

-0.0654 

0.0  35  1 

-0.8830 

-3.051  3 

0  .042  ? 

-0. 3  83  7 

-3 . 9  03  4 

0.0345 

-0.3355 

-0 . "  91 4 

3.0418 

0.1537 

.  O  1  'A  i 

0.0344 

0.2136 

3.9061 

9.0  41  R 

^  6971 

9.0472 

0.0350 

0.7634 

0.9797 

0.0  42  3 

1.2552 

2.0 974 

9. 0  35 7 

1 . 3  nQ7 

0.3741 

9.0  43  1 

1.7835 

9.1151 

9.0369 

1 .8708 

0.1137 

0.0445 

2. 3120 

9.193C 

3.0384 

2. 4095 

9.1  793 

9.9463 

2.  8746 

0.1719 

0.0  40  7 

2.9289 

0.1641 

3.0492 

3.4166 

0.2028 

9 . 043  1 

3.4  79  8 

9.2037 

0  .  C  5 1  2 

4 . 9  C  2  9 

0.2502 

3.0472 

4.053° 

0.2 3° 4 

C  .  054  6 

4. 529C 

0.2761 

9.  CSC  5 

4.  67Q1* 

9 . 2  5  Q  7 

0  .0  57  9 

5.0650 

9.3332 

0.0  54  6 

5. 1 330 

0.296  3 

9.0623 

5.6275 

0.3458 

9.0599 

5. 6894 

9.3366 

9.0679 

6.1656 

9.7706 

0.0651 

6.2  50° 

3 .2  "’6  0 

0 .0740 

6.7107 

0.4029 

0 .070  7 

6. 789° 

0.4912 

C  .  0  7°  2 

7.  2506 

9 .4249 

0.0765 

7.3267 

3.4336 

D  .086  3 

7.7930 

9.4571 

3.0842 

7 . 8  7°  1 

0.4563 

9  .  0  92  8 

8.  3412 

3.4871 

3.0919 

8.4152 

0.9823 

0  .  1  CO  2 

3.8710 

3.5143 

0 .0998 

8.  96,°7 

8.e  156 

9.1092 

o.  4  36  3 

0.5469 

0.  1093 

9.5042 

3.5454 

0.1189 

9.9823 

3.5750 

0.1187 

19  .  r>  337 

9 ,5638 

C .  126  5 

10. 5C87 

3.6002 

3  .  1 28  1 

10.6077 

3. 5 9° 2 

9. 1  399 

10.6730 

3.6111 

3.1315 

1C.  7465 

9.60  2  3 

0  .  1  4C  1 

MACH  NO. 
8  E/FT ( M) 

=1.026 
=  7.  99 

1  =  30 

MACH  NO. 
RE/  F^  (  M) 

=' , 699 
=  7 , 97 

1  =  38 

ALPHAS 

CL  B 

CD 

ALPHAB 

CL  B 

CD 

-4. 7  905 

-3.2935 

0.9 58 R 

-4. 8 1 3  7 

-0.2806 

6 .0587 

-4.  1  9°  1 

-0.2  57  1 

C.  0  537 

-4.21-3 

-6.2451 

0 .653  8 

-3.6  454 

-0.2139 

9.0495 

-3.6942 

-0.2232 

9 ,050  9 

-3 . 0994 

-0.1 868 

9 .046  7 

-3.1 32? 

-9.1 348 

0 . 0  47  7 

-2.5834 

-3.1 589 

0.0446 

-2.6095 

-3.1528 

9.0459 

-2.0126 

-9.1 266 

9.0424 

-2.0467 

-C  .  1 1  5  8 

0.0423 

-1 .4646 

-0.C91  3 

9 .C4C 7 

-1.5)27 

-0.0897 

C 

-0.921C 

-9.0527 

9 . 9  3°  5 

-0.9698 

-9.C524 

C.O  39  7 

-0.3598 

-9.9  218 

9.0390 

-0.4108 

-0.0177 

0.6392 

0.1860 

0.9170 

9.0390 

0.1 24 7 

0.0013 

0.03°  1 

0.7286 

0.8473 

C.C395 

6.6755 

0 . 0  36  3 

9.  0  34  6 

1. 2734 

0.08C9 

6.0  4C2 

1.2198 

0.0704 

0.04C  4 

1.8288 

0.H39 

0.0417 

1.7658 

0.1044 

0.0417 

2. 3594 

0.1356 

0.0  433 

2.3161 

0.1317 

0.0435 

2.9060 

8.1689 

0.0455 

2. 8603 

0.1664 

0.0458 

3.4693 

6.2154 

0 .048  7 

3.4195 

0.1993 

6 , 0  4  0  6 

4. 021 6 

0.2482 

0.C51 8 

3. 9578 

0.2260 

0.0  51  6 

4. 5726 

9.2812 

0.C559 

4.5072 

0.2652 

0  .0556 

5.0889 

0.3949 

9.0597 

5.061  6 

0.2  PR  8 

0.0604 

5.  6549 

0.3372 

0.0647 

5. 6048 

0.33C7 

0.0652 

6.  1990 

0. 3728 

0.07C3 

6. 1 409 

0.3559 

0.0  70  3 

6.7454 

0.4066 

0.0765 

6.7128 

0.3957 

0.0772 

7.2902 

6.4350 

0.0829 

7.2439 

0.4210 

0.0831 

7.8  328 

0.4617 

6.0897 

7.8092 

0.4529 

0.0906 

8.3771 

0.4909 

0.0973 

8. 3348 

0 . 4  75  8 

0.0975 

8. 9086 

0.50°6 

C.  1042 

8.8778 

0.4  99  5 

0.  1051 

9.4636 

0.5445 

C.  1 139 

9.4152 

0.5215 

0.  1129 

10.0169 

6.5734 

9.9788 

9.5548 

0.1232 

10.5501 

9.5968 

10. 5075 

0.5837 

0.  1331 

10. 7117 

9.6028 

C. 1349 

10.6944 

0.5945 

0.1369 
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TABLE  4 

AERODYNAMIC  COEFFICIENTS  OF  WING-BODY  MODEL  WB1  (Cont’d) 


MACH  NO.  =  ".956  M=3a 

P  E/FT (M ) =16.14 


ALPHAB 

CL  B 

CO 

-6.  2987 

-0. "41 7 

0  ,C8->7 

-5. 645C 

-" .  3C40 

"  .0467 

-5.0213 

-0.2593 

0 ."399 

-4.4442 

-0.2341 

0.C357 

-3 . 3794 

-0.1  O'M 

0 . C  31 3 

-3. 3123 

-0.1 846 

0.0287 

-2.7130 

-0 . 1 4  77 

C .0265 

-2.0709 

-C.101 5 

0.0230 

-1.5361 

-0.0896 

n  .(^21  " 

-0. 903  9 

-0.C54O 

(i/.irs 

-0.3069 

-0.0151 

0.C2r  4 

0.3096 

0.0257 

C  .0200 

0.8678 

0.C  505 

0.C2C7 

1.4642 

0.0322 

0.C21  5 

2.0  622 

C.  1  <"88 

C.0229 

2.6210 

'".1317 

0.0245 

3.2359 

",  1702 

0.0279 

3. 8334 

0.2Q^2 

0 .031  1 

4. 4232 

0.2295 

3. 0348 

5.0438 

0.2771 

9.04C4 

5.6326 

0.3  1C  5 

0.0459 

6.2771 

0.3555 

C  .053  6 

6.8431 

0.3846 

0.0  6C  4 

7.456C 

0.4184 

0.0683 

8. 0  5"  4 

0.4532 

C  .0772 

3. 6366 

9.4-»88 

0.0857 

9. 2270 

0.5  08  1 

0.0951 

9.7903 

0.5347 

0.  104  4 

10.3826 

0.5642 

0.1154 

10.8394 

0.5958 

0.1260 
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TABLE  5 

AERODYNAMIC  COEFFICIENTS  OF  WING-BODY  MODEL  WB2 


MAC  •!  10. 

=0.419 

'  < 

II 

MACH  94. 

VI  -  p  ^ 

FE/F?  (M) 

=  8.04 

0  5/  t>T>  f  \J) 

-  1  .  O  q 

ALPHAS 

CL  B 

CD 

ALPHAB 

CL  3 

Cl) 

-2. 8951 

-0.2806 

0.0197 

-3. "-72 

--.4173 

: . -  27  9 

-2.541'' 

-0.2414 

0.0166 

-2.6354 

-0.2747 

*  *  2 °  ? 

-2.2  419 

-P .2285 

0.01c  2 

-2.2319 

-  4 . 2  2  r  8 

**  O  'r  * 

-1 .9694 

-0.2138 

0.014  0 

-1  .  959'1 

-1.212? 

'  *  '  i  'i  r. 

.  •  • 

-1.6442 

-9.164  5 

0.4118 

-1.6434 

.  4  6  7 

-1  .2439 

-  0 .  ■*  2  O'  7 

0.0  1C  2 

-1  .  3196 

-1.14  6  2 

-  1 6  2 

-0. 9284 

-0.4  86  5 

0 .  CC°  3 

-1.04  19 

-0.1371 

••  ,r  14  7 

-0.6144 

-0.0  57  5 

0.0089 

-0. 6697 

-?. ; 754 

1  n  1  "3  . 

-0.2  9-'? 

-0.C330 

0.0086 

-0 . 1  99  5 

.-i/Of  q 

: .‘  0 1 4  4 

0.0368 

0  .  n  r>  2  4 

4.0485 

9.04IC 

9.71  ?  5 

.■'133 

n  3qtq 

0.4386 

6.0087 

0. 3556 

0.-457 

1  .  0  1  4  4 

0  ’.  6  996 

0.0  64  9 

O.C089 

3.6  75  5 

->  <■  7c  -> 

: .  ■- 1 2  ° 

1 .0  368 

0.  Ki  7 

0 .0096 

1.4  735 

4  .*  1  1  4  1 

. r  1 4  - 

1 . 3486 

0.1259 

0.010  4 

1.  3  51  1 

9.1  484 

-.•'754 

1. 7147 

0.1727 

C .C  120 

1 . 6  8  4  ( 

4.1774 

: .  c  1  {■  9 

2.0449 

0.2062 

0.0136 

2 .  ■'"'UP 

>  7  ,  =,  5 

.  r  1  O  5 

2.  3514 

0.2328 

4.0152 

2.3436 

7  !  4  4  P  8 

6.0214 

2.6960 

4.2738 

0.0177 

2.6734 

4.2784 

j ,  2  2  3  8 

3.0755 

4.3269 

0.0  216 

7  "■''K 

4.3  145 

4 . C  27  0 

3.3997 

0.3539 

C .024  ? 

3*.  3  64  5 

1.3689 

■j  .  ^  4 1  "> 

3. 6921 

4.3442 

0.0268 

3.  6  76-1 

4.7  846 

0.0 346 

4.0467 

0.4171 

C.031 5 

4.  002  5 

9.4122 

4.0382 

4.  9705 

4.4518 

0.0359 

4.303? 

4.4376 

0.0421 

4.6778 

4.4753 

0.0398 

4.6337 

0.4695 

0.0466 

5 . 0  0 1  9 

4.54^5 

4.0449 

4.  9774 

0.5  06  2 

0 .0  522 

5. 3080 

0. 531 9 

0 , 4  4Q  2 

5.2764 

9. 5260 

9.0  56  3 

5.  6  39  6 

0.5644 

0.C550 

5.  591  6 

0.5556 

0.0617 

5.963  3 

4.5928 

0.0605 

5. 8656 

0.5639 

0.0654 

6.0917 

4.5  02  4 

0.0620 

*  ACI1  NO. 

=0.478 

1  =  29 

MACH  NO. 

=  9.  979 

OO 

CM 

II 

*-• 

RE/FT  (*) 

=  "'  .  04 

RE/FT(M1 

=  9.44 

ALPHAB 

CL  B 

CD 

ALPHAS 

CL  3 

CD 

-7,0434 

-  7.3  60  8 

0.0  27  2 

-2. 9481 

-0.2704 

0 .029  5 

-2.724? 

-9.3219 

C  .023  9 

-2.  6123 

-0.23Q9 

0.0270 

-2.3466 

-0.2  72  6 

0.019  3 

-2. 3042 

-4.2224 

0.0254 

-1. 9492 

-0.2086 

0.0156 

-1.9852 

-0.1 939 

0.0734 

-1.6  93  6 

-4.1 879 

0.0141 

-1.6799 

-0.1613 

4.021 7 

-1 . 3451 

-2.1 554 

0.0124 

-1.3774 

-0.1442 

0.02-3  5 

-1.-146 

-0.1216 

4,011" 

-1. 0019 

-0.0956 

0.C190 

-9.6  57° 

_r  "  7  «  1 

•  . 

'  0  "  4  C 

-0.6663 

-O.O519 

0. 0181 

-2  .  3  36? 

-0  .'-4  6  1 

9  !  c  0  4 

-0.3  52  0 

-9. 0328 

0.0177 

0 . 0  4  ) 0 

7  -.r  or 

44  0  0  r>  ^ 

•  •  -  -  - 

-0.0459 

-0.0069 

0.0176 

0.  3933 

7  *  6  (I  8  6 

c .  rc°4 

0.2677 

0.0171 

0.0178 

0.6  931 

? .4766 

i.PV 

0.6078 

0.0559 

0.0181 

1 . 4 345 

4.1 3P  7 

0 .  c  1  i  2 

0.9300 

0.0844 

0.0187 

1.3933 

0.1 566 

:  .4125 

1. 2736 

4.1208 

0 .0199 

1. 7380 

•9.  1°6  8 

9.0146 

1.6395 

9.1741 

0.021 9 

2. 0756 

2. 4266 

9.2 3? 5 

n  ^  n 

".0167 
°  .0  19  6 

1.9537 

2.2666 

0.1937 
0.221  1 

2*.  0^5  0 

2.6997 

41284? 

0 .021 2 

2.  6106 

0.2587 

0.0  27  3 

3.  0  57° 

4.? 35 6 

0.0251 

2.8971 

0.2798 

C.0291 

3.3711 

0.3592 

0.0279 

3. 2288 

0.3107 

0.0  319 

3.  6462 

4 . 37” 2 

.  0  3C  1 

3.  5754 

0.3496 

0.0356 

4 . C  592 

9.4380 

'7.036  8 

3.8902 

0.374Q 

0.0387 

4. 3782 

4.4676 

r  .C4C  9 

4.1904 

0.3965 

0.0420 

4. 701 4 

4.4476 

C.CU53 

4.5819 

C.4514 

0.0480 

5.0  1°2 

9.524  2 

C  .  0  5C  3 

4.8364 

0.4588 

0 .050  0 

5.3396 

9.5559 

0.0558 

5.2066 

0.5055 

0.0565 

5. 6241 

0.5666 

0.0  59  3 

5.5023 

0.5254 

0.0607 

5.9894 

0.6107 

4.0671 

5.8666 

0.5705 

0.0677 

6.1159 

4.6140 

0.4689 

I 
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TABLE  5 

AERODYNAMIC  COEFFICIENTS  OF  WING-BODY  MODEL  WB2  (Cont’d) 


MACH  NO. 
PE/FT  ( 1) 

=  o. 80  i 
=  9.92 

1=29 

MACH  NO. 
RE/ FT  (  M) 

=P.874 
=  7.  99 

1  =  28 

ALPHA B 

CL  3 

CD 

ALPHAB 

CL  B 

CD 

-2. 7220 

_/>  i  oo  ■> 

C .013  3 

-2. 777t 

-0.2481 

0.0152 

-2 . 4 1 1 C 

-0 !l757 

0.0120 

-2. 4838 

-0.2371 

0.0143 

-2.0509 

-9.1456 

0 .0109 

-2.0  95? 

-0.1854 

0.0H9 

-1.9  322 

-9.1425 

0.0108 

-1.7659 

-0.1462 

0.0106 

-1.9134 

-9 . 1 1 1  0 

^ . 0997 

-1 .4997 

-0.1281 

9.01C2 

-1. 221 5 

-0. 1 r95 

9.0093 

-1.1157 

-0.0844 

9.0090 

-P.  8  99 ^ 

-0  .''690 

0.0986 

-0.8  56  8 

-0.0772 

0.0088 

-9.6309 

-0.0523 

0.9093 

-0.599'’ 

-0.9361 

0.0092 

-0.3139 

-° . C  28  4 

0.0080 

-0.2036 

-0.0169 

0.0C81 

0.0059 

-0.9018 

0.0089 

0.1278 

0.0174 

C .0980 

0.2794 

9.C139 

9 .0080 

0.4380 

0.0395 

0.0093 

0.6122 

9.0486 

9.0084 

0. 7328 

0.P621 

0.0087 

n  ,  q  04  6 

o . 0  692 

P  .008  9 

1.0143 

0.9737 

0.0090 

1.2448 

9. 0986 

0.9994 

1.3471 

9. 1056 

0 . 0097 

1. 5685 

0.1398 

0.010  2 

1.6  72  3 

0.1  33  6 

9.  Cl? 2 

1 .8532 

9.1487 

0.0  1C  5 

1 . 9676 

0.1521 

0.0107 

2. * 657 

C.175C 

p.0117 

2.3537 

0.2100 

0.0130 

2.4401 

0.1878 

0.0124 

2.6  53  4 

9.2296 

0.0143 

2.  7642 

9.2101 

9.0139 

2.9764 

0.2606 

0.0162 

3.9664 

C.2435 

0.0155 

3.  2988 

0.2894 

0.0183 

3.  3-766 

3.2656 

0.0174 

3.  6  310 

9.3260 

0.021 2 

3.6570 

0.2828 

0.0189 

3.9642 

9.3546 

0.0239 

4.0959 

0.3210 

9.0  222 

4.2632 

0.3744 

9.0268 

4.2984 

9.3366 

3.0  24  4 

4. 5873 

O.405O 

P  .  0  30  9 

4.6144 

0.3699 

0.0280 

4 .  8  54  4 

0.4127 

0.0334 

4. 8972 

9.3790 

9.0  30  0 

5. 1 653 

0.4353 

0.0369 

5.  1640 

0. 38-74 

0  .  C  32  4 

5.4405 

0.4455 

0.0394 

5.477C 

0.4286 

9  .  C  37  8 

5. 7396 

0.4793 

0 . 0  44  7 

MACH  NO. 

=0.849 

1=29 

MACH  NO. 

=0.999 

1=29 

PE/FT (M) 

=  8.03 

RE/  F'1’  { M) 

=  8.01 

ALPHAB 

CL  B 

CD 

ALPHAB 

CL  B 

CD 

-2. 704  0 

-9.2375 

0.0152 

-2. 8991 

-0.2744 

0.0181 

-2. 3363 

-0.1899 

0.0124 

-2.5154 

-0.2183 

0.0145 

-1.9861 

-P.15R5 

C  .  0  1 1  5 

-2.1 785 

-0.1981 

9.0128 

-1 . 7643 

-0.1543 

9.0109 

-1.8615 

-0.1674 

0.0115 

-1.4246 

-9.1179 

P.CP98 

-1 . 5689 

-9.1358 

0.C1P5 

-1.097Q 

-9.9889 

0 .  p  C  8  9 

-1.2745 

-0.1292 

0.0097 

-0.7502 

-0.9 UUQ 

p  .0080 

-O.9317 

-0.0870 

0 .0999 

-0. 5022 

-0.0442 

0.0081 

-0.6224 

-0.0563 

0.0005 

-0.1 549 

-3.0126 

p.0978 

-0.2906 

-0.9258 

0.CC92 

9.1237 

9.P098 

0.PP78 

0.0611 

0.0103 

0.0081 

0.4335 

0 . P  28 1 

9.0980 

0.3578 

0.0349 

0.0082 

0 . 7530 

P.C6P7 

9.0984 

0.6900 

0.0649 

0.0086 

1.0772 

0.0845 

n  .onoo 

9.  9967 

9.0939 

0.0090 

1.  7  67  6 

0.1030 

0.P096 

1 . 2855 

0.1020 

0.0095 

1.6892 

0.1  369 

0.019  3 

1.6465 

9.1 495 

0 . 0  IP  8 

2.  p  129 

P.1644 

0.0111 

1 .9698 

P.1  834 

0.0120 

2. 3093 

9.1 857 

0.0122 

2.3044 

9.2142 

0.0134 

2. 6985 

0.208 2 

0.0134 

2.6778 

p.2659 

0.0159 

2.9436 

3.2438 

0.0151 

2.9659 

0.2830 

0.0176 

3. 2281 

9.2512 

0.0  164 

3.  26C7 

0.3C12 

0.P193 

3.5711 

3.2929 

C.C1P6 

3.5896 

9.3298 

9.0220 

3.  854° 

9.3093 

0.0206 

3.9470 

0.3776 

0.0263 

4.2106 

9,351 3 

0 .024  4 

4.2565 

0.4026 

0.0295 

4 . 5 2 5n 

p.3720 

0.0270 

4.  5637 

0.4222 

0.0328 

4.8442 

3.4039 

0.07^7 

4. 8670 

9.4447 

0.0366 

5.1197 

0.4125 

0 .0330 

C . 0379 

5.  1 91 4 

0.4723 

0.041  1 

5.4253 

0.4367 

5.5022 

0.4962 

0.C456 

5. 6431 

P.4404 

0.0390 

5.  7810 

9.5049 

0 .0489 

5.9127 

0.5C37 

0.0502 
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TABLE  5 

AERODYNAMIC  COEFFICIENTS  OF  WING-BODY  MODEL  WB2  (Cont’d) 


HACK  MO. 
RE/FT  (M) 

=  0.  998 
=  7 ,  qq 

*=28 

MACH  MO. 
RE/FT  (M) 

=1.050 

=  7.46 

!i=2  8 

ALPHAB 

CLB 

CD 

ALPRAB 

CLB 

CD 

-2.7925 

-0.2530 

0.0293 

-2.8654 

-0.2631 

C.034  1 

-2.3897 

-0.2044 

0.0260 

-2.4941 

-0.2217 

0.031  1 

-2.1122 

-0.1  967 

0.0250 

-2. 1 665 

-0.1 899 

0.0289 

-1.8252 

-0.1772 

0.0238 

-1.8327 

-0.1 554 

0.0270 

-1. 5064 

-0.1384 

0 .022  0 

-1.5649 

-0.1 392 

C . 0260 

-1.1508 

-0.1030 

0.02C  7 

-1.1715 

-0.0936 

0.n243 

-0.8326 

-0.0774 

4.0200 

-0.9047 

-o.oe27 

C  .023  9 

-0.5540 

-0.0 56 3 

0 .019  5 

-0.6 081 

-O.C583 

0.C234 

-0. 1641 

-0.0053 

0.014C 

-0.2715 

-0.0274 

0.023  1 

0.1395 

0.0165 

0.0191 

0.0267 

-0.0078 

0.0231 

o.  4  72  6 

0.CU48 

C.C19  3 

0.4012 

0.0362 

C  .r23  2 

0. 7489 

0.065  8 

0.0194 

4.7164 

0.4682 

0.0236 

1.0016 

0.0460 

0  .  C  20  5 

1.0652 

0.1CC2 

0.0243 

1 .4029 

0.1242 

0 . 021 4 

1.  3810 

0.1 278 

0.0253 

1.7240 

C.15C6 

0.0227 

1 .686  5 

3.1 499 

0.0266 

2.0  521 

0.1  838 

0.0244 

1.9845 

0.17^2 

0.C278 

2. 3890 

0.2162 

0.0260 

2.  3335 

0.2062 

0.0297 

2.7154 

0.2520 

0 .028  3 

2. 6715 

0.2449 

0.0321 

3.0569 

0.2848 

4.0  310 

2.9863 

n .  4  6  4  o 

0.0341 

3.3908 

0.3147 

C.0340 

3. 3295 

0.3084 

0.0377 

3.6967 

0.3412 

0.0  366 

3.6379 

0.3300 

0.0440 

3.9683 

0.3558 

0.C389 

3. 9409 

0.3508 

C.C42  6 

4.  3553 

C.4085 

C.C44  3 

4.2854 

0.3856 

0.046  5 

4.6437 

0.4220 

0.0470 

4. 5949 

0.4114 

0 . 0  5C  1 

4. 9568 

0  .4475 

C .0510 

4.8978 

0.4324 

C.0536 

5.3213 

0.4965 

0.0573 

5.2034 

4.4522 

0  .  C  57  4 

5.6197 

0.5146 

0.0616 

5.5497 

0.4  92  9 

0.0632 

5.9651 

0.5535 

0.0675 

5. 8688 

0.5142 

0.C682 

MACH  NO. 
RE/FT  (M) 

=1.025 
=  8.02 

*=28 

MACH  NO. 
RE/FT  (*) 

=  ‘'.099 
=  7.  °° 

M  =  2  8 

ALPHA  B 

CLB 

CD 

ALPHAB 

CLB 

CD 

-3.  0360 

-0.2707 

0.0332 

-2.9086 

-0.2515 

0.C328 

-2.7145 

-0.2532 

0.C31  3 

-2.6985 

-0.2465 

0 .C  31 8 

-2.3  560 

-0.2075 

0.0283 

-2.  2520 

-0.1832 

0.0277 

-2.0529 

-0.1  86  9 

0.0270 

-2.0000 

-0.1 679 

0  .  °26  5 

-1.7  54  8 

-0.1558 

0 .026  3 

-1 . 7252 

-0.1487 

0 .0263 

-1.4215 

-0.1337 

*.0240 

-1.3312 

-0.1193 

0.0239 

-1.1C59 

-0.1047 

0.023C 

-1.0241 

-0.081 0 

0.0226 

-0.7217 

-0.0472 

0 . 0  21 7 

-0.7  303 

-0.0595 

0.0220 

-0.4318 

-0.0320 

0 . 0  2 1  5 

-0.4155 

-0.0373 

0.021 7 

-0.1565 

-0.0244 

0.021 6 

-0.0922 

-0.0063 

0.021 u 

0. 1912 

0.0172 

0.021 5 

0.2280 

0.0199 

0.021  6 

0.4933 

0.0408 

0.0217 

0. 5476 

0.0472 

0.0218 

0.8609 

0.0822 

0.0223 

0.8512 

o.07i ° 

0.0223 

1. 1689 

0.1C48 

0.0230 

1. 1 540 

0 ,°905 

0.0230 

1.4703 

0.1 31 2 

0.0240 

1 .4880 

0.1210 

C.C241 

1.8088 

0. 1602 

0.0254 

1.8284 

0.1657 

0 .0256 

2.  1231 

0.1959 

0.0268 

2.1708 

0.1924 

0.0276 

2.4473 

0.2194 

0.0287 

2.485° 

0.2176 

0 . 0  2°  2 

2.7780 

0.2489 

0.C3C8 

2.7831 

0.2350 

0.0308 

3.  1091 

0.2 84 Q 

0.0  33  7 
0.0372 

3. 1 184 

0.2677 

0.C  335 

3.4569 

0.3237 

3.4534 

0.3017 

0.0  36  7 

3.7992 

0.3  60  2 

0.0406 

3.767° 

0.3246 

C.0394 

4.0944 

0.3786 

0.0434 

4. 0864 

0.3541 

0.043  1 

4.4568 

0.4253 

0.0486 

4.  459° 

0.3988 

0.C48  1 

4. 7761 

0.4501 

C  .0522 

4.7396 

0.4130 

C.0509 

5.0975 

0.4782 

0.0565 

5.0614 

0.4  40  3 

0.0551 

5.3961 

0.4960 

0.0602 

5.3990 

0.4752 

0 .060  3 

5.  7332 

0.5383 

0.0666 

5.  6591 

0.4879 

0.F635 

FIG.  3:  CROSS-SECTIONAL  AREA  DISTRIBUTIONS  OF  THE  MODELS  WITH  AND  WITHOUT  LIFT 
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FIG.  5:  VARIATION  OF  CL  WITH  a  AND  Cp  WITH  CL  FOR  THE  WING-BODY  MODEL  WB1 

AT  DIFFERENT  MACH  NUMBERS 
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FIG.  6:  VARIATION  OF  CL  WITH  d  AND  CD  WITH  CL  FOR  THE  WING-BODY  MODEL  WB2 

AT  DIFFERENT  MACH  NUMBERS 


WING  BODY  W8I 


B2,  a  =-0.05° 


FIG.  9:  SCHLIEREN  FLOW  VISUALIZATION  FOR  WING-BODIES  WB1,  WB2  AND  BODY  OF 
REVOLUTION  B1  AT  ZERO  ANGLE  OF  ATTACK  (NOMINAL),  M„  =  0.98 
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FIG.  10:  STATIC  PRESSURE  DISTRIBUTIONS  NEAR  THE  UPPER  AND  THE  LOWER  WALLS 
OF  THE  TEST  SECTION  FOR  THE  WING-BODIES  WB1,  WB2  AND 
THE  BODY  OF  REVOLUTION  B1  AT  ZERO  ANGLE  OF  ATTACK,  M„.  =  0.98 


*T  ZERO  LIFT  WITH  TRANSONIC 
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FIG.  12(b):  SURFACE  FLOW  VISUALIZATION  FOR  THE  WING-BODIES  WB1  AND  WB2 

AT  ANGLE  OF  ATTACK,  M„,  =  0.96 


MACH  NO.  =0.8 


FIG.  14:  VARIATION  OF  CDCL_  WITH  C  2  FOR  THE  WING-BODY  WB1  AT  DIFFERENT  MACH  NUMBERS 


FIG.  15:  VARIATION  OFCnC,  _  WITH  C,  2  FOR  THE  WING-BODY  WB2  AT  DIFFERENT  MACH  NUMBER 


VARIATION  OF  THE  INDUCED  DRAG  FACTOR  k  WITH  MACH  NUMBER  FOR  THE  WING-BODIES  WB1  AND  WB2 


ly  +  l)M„T-b 


FIG.  18:  DRAG-RISE  FOR  THE  WING-BODIES  WB1,  WB2  AND  THE  BODY  OF 
REVOLUTION  B2  AT  THE  DESIGN  LIFT  CONDITION 


FIG.  19:  SURFACE  FLOW  VISUALIZATION  FOR  THE  BODY  OF  REVOLUTION  B2 

AT  MACH  NUMBER  0.96 


STATIC  PRESSURE  DISTRIBUTIONS  NEAR  THE  UPPER  AND  THE  LOWER  WALLS  OF  THE  TEST  SECTION 
FOR  THE  WING-BODIES  WB1  AND  WB2  AT  THE  SAME  LIFT  PARAMETER,  M„  =  0.94 


APPENDIX  A 


SIMILARITY  PARAMETERS  OF  THREE-DIMENSIONAL  TRANSONIC  FLOWS 


Similarity  parameters  of  a  physical  problem  can  be  derived  from  the  governing  differential 
equations  and  the  boundary  conditions.  For  three-dimensional  transonic  flow  with  small  disturbances 
the  parameters  come  about  naturally  in  the  derivation  of  the  small  disturbance  equations.  With  the 
flow  taken  to  be  uniform  upstream  in  the  x-direction ,  the  velocity  components  and  the  density  of  the 
flow  can  be  written  for  a  small  disturbance 

u*  =  U_  (1  +  Au) 
v*  =  U„  Av 

.  (Al) 

w  =  U_  Aw 

P *  =  (1  +  Ap) 

where  A  is  a  perturbation  parameter  to  be  determined  from  the  small  disturbance  condition,  i.e.  the 
boundary  condition  on  the  body .  The  co-ordinates  are  scaled  as  follows 

x*  =  Cx 

y*  =  by  (A2) 

z*  =  c  z 

where  C,  b  and  c  are  the  typical  lengths  of  the  body  in  their  respective  co-ordinate  directions.  If  the 
body  surface  is  given  by  the  expression 

z*  =  5  f*  (x*,  y* )  (A3) 

where  6  is  small,  so  that  the  requirement  for  small  disturbance  is  satisfied,  we  can  show  that 

A  ~  0(6)  (A4) 

The  governing  equations  chosen  for  the  derivation  of  the  small  disturbance  equation  are 
the  conservation  equations  of  mass  and  energy.  These  equations  satisfy  automatically  the  discon¬ 
tinuity  conditions  of  the  shock  wave  which  may  occur  in  the  flow  field. 

(P*  u*)  ,  +  (p*  v*)  ,  +  (p*  w*)  ,  =0 

x  y  z 

,2  a*2  (A5) 

—  (u* 2  +  v* 2  +  w* 2 )  +  -  =  —  U* 2  + - 

2  7 -  1  2  7-1 


With  the  perturbation  velocities  Equation  (Al)  and  eliminating  the  density  with  the  energy  equation 

.2 

(A6) 


P*  ,  r  1  1  1  r  1  m2 

—  =  !  +  Ap  =  1  -  M2  |Au  +—  (Av)2  +  -  (Aw)2J  +  [(2  -  7)  M2  -  lj  (Au)2 


the  perturbation  equation  for  the  velocity  is  derived  up  to  the,  nrrf« 


KaCSDlNfc  PJU OS  BUMMWT  FILMED 


T 


A  ,  A  A  M2  a2  ,  ,  A2  M2  , 

-  (1-  M2)UX  +  -  vy  +— Wz  -  7  (v2  +  w2  )x  -  —  (7+1)  —  (U  )x 


M' 


A2(—  u  v  +—  u,  w)  -  M2  A2  f  —  uv„  +  —  uw  )  =  0 
\b  y  c  z  /  ~  \b  y  c  V 


Assuming  that  the  flow  is  irrotational  and  introducing  a  velocity  potential 
<i>*  =  U_  x*  +  0*  =  U„C  (x  +  A0) 


U 


with 


90  C  90  2  9 0 

u  =  —  ,  v  =  77  ,  w  =  -  — 

9x  b  9y  c  9z 

The  velocity  equation  can  be  written  in  terms  of  0  as 

[(I-  MJ  "  A(7+  1)M„0x]  0xx  +  p  ^yy  +jr  tzz  =  <Py  +  ^]x 

Tc2  c2  1 

+  AM2„  [p-  0yy  +  p-  0zzj0x 


For  a  configuration  with  small  or  moderate  aspect  ratio,  the  length  scale  c  is  taken  as  b,  and  Equa¬ 
tion  (A9)  can  be  rearranged  into  the  following  form 


1-  Ml 


-  0V 


(7  +  1  )M2  A  rx 


C2 


(7  +  1)M„  A  b2  yy 


7  +  1  b2 


— T  (0  +  0  )  = - —  — r  (0  +  0  ) 

*2  a  k2  '+yy  ^zz'  ^  a.  i  k2  yvy  vz’ 


i  e2 

+  y  +  1  b2^yy  +  ^ZZ^X 


(a: 


The  perturbation  parameter  now  has  the  form 


A  = 


=  rX 


(a: 


6b  b 

where  r  =  —  and  X  =  —  are  the  thickness  and  the  swept  back  parameters  respectively.  The 
boundary  condition  is  therefore 


90  _  9f(x,  y)  a 
9z  9x  r 


y  =  0 


(A] 


With  the  perturbation  equation,  Equation  (A10),  and  the  boundary  condition,  Equation  (A12),  the 
similarity  parameters  can  be  identified  as 


1-  Mj 

(7  +  1)  Ml  rA  ’ 
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a 


(7  + 1)M  IrX3,  (7+l)X2'  — 


These  parameters  appear  as  coefficients  in  the  equations  defining  the  order  of  magnitude  of  their 
related  terms.  They  can  be  identified  as  follows:  the  first  one  is  the  transonic  similarity  parameter  K, 
the  second  one  the  square  of  the  ratio  of  the  length  scales  for  the  outer  and  inner  regions,  e,  the  third 
one  the  swept  back  parameter,  related  to  T,  and  the  fourth  one  the  normalized  angle  of  attack, 
related  to  a ,  (see  Eqs.  (1)  and  (4)  in  Section  2).  Because  of  the  retention  of  the  right-hand  side  in 
Equation  (A10),  these  parameters  are  not  the  same  as  those  of  the  classical  similarity.  Different  forms 
of  small  disturbance  equations  have  been  used  in  the  literature  on  three-dimensional  flows*1 12 .  They 
are  basically  the  same  as  Equation  ( A10)  except  that  some  terms  in  the  right-hand  side  are  dropped 
according  to  the  author’s  preferences.  If  all  terms  at  the  right-hand  side  of  the  equation  are  dropped, 
the  classical  three-dimensional  equation  is  obtained  and  the  classical  similarity  law  follows*102. 
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APPENDIX  B 


EVALUATIONS  OF  FUNCTIONS  T(x),  E(x)  AND  F(x)  FOR  A  DELTA  PLANFORM 


For  a  planar  wing  at  angle  of  attack  the  cross-flow  potential  given  by  the  slender  wing 
theory  is(12) 


V  =  (a2  -  y2)1/2 

where  a(x)  is  the  contour  of  the  leading  edge. 


(A13) 


The  functions  T(x),  Equation  (6)  in  Section  2,  can  be  evaluated  with  the  potential  given 
in  Equation  (A13) 


T(x)  f  f  Mx,s)]x  Mx,y)]x  fin  |  — ^ —  |  ds  dy 

47T  -°°  y  s 


a2  a  2  .a 


*  -2  _  “2  vi/2 


7T  .a  (az  -  y2 ) 
The  integral  with  respect  to  s  can  be  evaluated  as 


2  -  «2vl/2 


(a-8  -  s') 


fin 


ds 


T  =  ra  _ l _  fin _ I _ ds 

1  -a  (a2  -  S2)1^2  I  y  "  s  | 


.  .  l  fin  |  y  -  s  , 

(fin  a)  7r  -  /  — — ds 


J-1  (1  --s2  )1/2 
(fin  a)  7 r  -  I2 


where  y  =  —  and  s  =  — 
a  a 


To  evaluate  I2  we  differentiate  I2  with  respect  to  y 


dl 


2  1  1 


=  / 


dy  ''-I  (1  -  s2)1/2  y  -  s 

Let  s'  =  cost?,  y  =  cosfij  then 


ds 


^  *2  n  d  (? 

=  /  — - -  -  0 


dy  o  cost?  j  -  cos# 


Hence  I2  =  constant. 


At  y  =  0  the  integral  I,  reduces  to 


_  fi  gnjaj 

A  9  J  “  __n "  i  /->  uS 


12  '-1  (1  -¥2),/2 

Since  the  integrand  is  an  even  function  of  J, 


[y-  s 


dy  (A14) 


ppy-gnitA  p^OS  BLAtK-NOV  TllMSD 


if  kn  2 


I  k'n  I  s  | 

U-  ■  2  {,  ,1  ,V-'  “* 


Thus 


And 


1,  =  n  t!n 


T(x) 


2  2 

a*  a 1  a 


(t) 

,tn(l) 

J  '  2  _  „2,l/2 


4n  (a  -  y  ) 

/ 

w  8n 

The  integral  E(x),  Equation  (7)  in  Section  2,  is  evaluated  as 


dy 


'(f)a2a'! 


E(x)  =  —  /"/  Wx,s)]  Mx,y)]y  £n 
47T  —  y  * 


-V  y 


IT  (a2  -  y2  )1/2 
The  integral  with  respect  to  s  is  evaluated 


/* 


(a2  -  s2 ) i/2 


y-  s 


£n 


ds  dy 


y-  s 


ds 


dy 


>.  "  C  *”  —L~  ds 


I  y-  s  | 


=  -  a 


l  s 

/  - — Cn  I  y  -  s  I  ds 
-i  (1  -  ss)l/* 


•  -  a  I, 


To  evaluate  the  integral  I2  we  again  differentiate  I2  with  respect  to  y 

1 


dii  = 


dy  U  (1  -  s2)1/2  |  y  -  s 


ds 


Substituting  s  =  cos#,  y  =  cos#, 


Thus 


^  *2  fO  cos# 
dy  a  I  cos# ,  -  cos#  | 

=  -7T 

I,  =  -jry 


(A15) 


(A16) 
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and 

I,  =  Try 

The  integral  E(x)  becomes 

E<X)  ‘  r  C  (.J-y2)i7J  <’ry»d>' 

=  I  2 

2  3  (A17) 

Finally,  the  integral  F(x)  is  evaluated  as 

FM  =  r  Mx,y)]  dy 

-oo 

=  2  /  (a2  -  y2 )  dy 

-a 

=  (A18) 

i 

; 

i 


-57  - 


APPENDIX  C 


BOUNDARY  CONDITION  AT  THE  TUNNEL  WALL 


The  linearized  boundary  conditions  for  flow  past  a  perforated  wall  is  given  as(19) 


*  1  .  A 

u  ±  —  v  =0,  y=±h 

P* 


(A19) 


where  P  is  the  porosity  factor  depending  on  the  physical  configuration  of  the  wall.  The  plus  and  the 
minus  signs  in  the  condition  refer  to  the  upper  and  the  lower  wall  respectively.  For  the  side  walls,  the 
velocity  component  w  will  replace  v  in  the  equation  and  the  treatment  is  identical  to  the  other  walls. 
Introducing  the  velocity  potential,  Equation  (A8),  the  condition  can  be  written  as 


30  j__e_  30  _ 

3x  p*  b  3y 


(A20) 


In  terms  of  the  similarity  variables  defined  in  Equation  (1),  we  have 


30  1  30 

—  ± - -  =  0 

3x  p  dr) 


(A21) 


where 


Se  "  [(7  +  DMjr  X]"2 


(A22) 


Thus  for  flow  similarity  at  the  wall  it  is  necessary  to  have  P  the  same  for  both  models.  Since  7  X  is  the 
same  for  the  models  according  to  the  first  condition  of  Equation  (26),  the  scaled  porosity  coefficient 
P  will  have  the  same  value  if  the  tests  are  performed  at  the  same  Mach  number. 
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APPENDIX  D 

ESTIMATION  OF  OFF-SETTING  FOR  CL  vs  a  AND  CD  vs  CL  DATA 


The  off-setting  of  the  CD  vs  CL  polar  from  the  zero  lift  axis  is  due  to  the  asymmetry  of 

the  model  about  the  central  plane.  With  the  angle  of  attack  range  of  -  3°  to  3° ,  the  data  are  fitted  with 
a  quadratic  in  the  least  square  sense, 


CD  =  A  +  BCL  +  CCL2 


The  minimum  of  CD  is  thus 


and  the  corresponding  CL  is 


=  A 


B2 

4C 


2C 


(A23) 


(A24) 


(A25) 


The  angle  of  attack  shift  is  taken  to  be  the  value  of  a  at  CL  by  a  linear  fit  through  the  data  of 

m 

CL  vs  a. 


CL  =  D  +  Ea 


Thus 


a 


The  data  are  then  presented  in  the  following  forms 


(A26) 


(A27) 


a  =  a  -  a 

m 


(A28) 


PRKSDlwO.  PMOS  BUMC-ttOT  71UCP 


